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ABSTRACT 


This study Bro lenes a theoretical investigation of 
the time-dependent, coupled conduction-convection heat 
transfer problem of a heated radial fin rotating simul- 
taneously with its surroundings at a constant angular 
velocity about a horizontal axis. Mathematical formulation 
reveals that this is a periodic problem governed by a set 
of unsteady-state, nonlinear, partial differential equations 
which have neither initial conditions nor synchronous 
solutions. 

Through an order-of-magnitude analysis, six para- 
meters appear in the dimensionless governing equations and 
boundary conditions. The Prandtl number Pr is an indication 
of the relative importance of viscous dissipation and heat 
COROUCtiOn in+ the. fluidsfiow;4thel Froude ntumber®Fr-is»the 
ratio of the centrifugal ee ceiereedon to the gravitational 
acceleration; the geometric parameter y is the ratio of the 
Sune onothie lOstheedrsvanee between theViin tip and thetcenter 
and K, are the 


is 2 


measures of the transient responses in the fluid and solid, 


of rotation; the response parameters K 


respectively; the sixth parameter C is an indication of the 


degree of coupling between the conduction and free convection 


heat transfer problems. 
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vViate the ‘effect of Prandtl number on the solutions. ‘The 
Prandtl number of air, i.e. Pr=0.72, is then used through- 
out the analysis. 

The fluctuating amplitude of the periodic solutions 
decreases inversely with the Froude number. [In the limiting 
case of very large Froude number, the guasi-steady-state 
solutions, which are independent of the parameters Ky and Kor 
are obtained. 

These quasi-steady-state solutions are presented over 
tee@range sot «0.1<ys0.5 and 0.01<C<10..' The results are in 
good agreement with those of previous work. 

Perturbation solutions are used to study the periodic 
heat transfer phenomena of small fluctuations. The investi- 
gation is conducted over the range of 10<Fr<100, 0.01<C<l10, 
OWw<y-<0: 4570). ToK 


<k0'O' > Fand 0.1<K,<1000. In general, periodic 


i 
fluctuation in the velocity solutions is appreciable when 
Fr=10, but negligible when Fr=100. On tie other hand, no 
significant fluctuation is present in the temperature solu- 
Evons $ 

The periodic solutions under small and moderate fluc- 
tuations are obtained by directly integrating the unsteady- 


state governing equations for this heat transfer problem. 


Because of the expensive numerical work required in the 
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approach of direct integration, solutions are only presented 
Over the range of 2<Fr<10> FOr Y=O. pyro 0. OL Kiel, and 
K,=100. The fluctuating amplitude of the solutions decreases 
inversely with the Froude number (Fr); however, the total 
heat transfer rate increases with the latter. 

Analysis also shows that the cooling effect of a 


1/4 times that of a conven- 


BOcating, titivis approximately. Fr 
tional stationary fin. In other words, one can double the 
SlLLErGlency Olea estatzonaly Cooling fin ti tt as. technically 
and economically permissible to rotate the system to l6ég 
(weet Fr=16)) 

Study of the solution applicability reveals that if 
the critical Rayleigh number for a stable laminar free con- 
vection flow is aye the solutions presented in this work 
are applicable to a rotating system with the centrifugal 
acceleration NO Gfeacer than 500%, <Showildy thesactural 
critical Rayleigh number obtained from the stability study 
of this periodic heat transfer problem be greater than Vel 
the presented solutions may have some application to the 
high speed turbines. 

A rational approach has been developed to solve this 
periodic heat transfer problem governed by the diffusion- 


type, nonlinear partial differential equations, of which 


the initial conditions are missing and the synchronous 
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solutions under the finite-amplitude fluctuations are not 
obtainable. This approach can be used as a guideline to 
solve other periodic problems, which occur in the physical 
world but remain untouched because of the large fluctuations 
in the non-synchronous solutions or because of the missing 


invtial conca tions. 
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Fin Effectiveness of Quasi-Steady Solutions 
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the Boundary Layer Problem 
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Effect of y on the Perturbed Boundary Layer 
Solutions 
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Perturbed Boundary Layer Solutions for 
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Perturbed Boundary Layer Solutions for 
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First-Order Perturbation Temperatures Along 
Solid Fan 


Periodic Variation in Perturbed Heat Transfer 
Rate 


Periodic Variation in Nusselt Number of 
Perturbation Result 
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Periodic Boundary Layer Solutions Obtained 
from Perturbation Approximation for Fr=10 


Periodic Boundary Layer Solutions Obtained 
by the Direct Integration Method for Fr=10 


Periodic Boundary Layer Solutions Obtained 
by the Direct Integration Method for Fr=2 
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Periodic Fin Temperature Obtained by the 
Direct Integration Method for Fr=2 


Periodic Variation tn the Total Heat 
Transrer Rate of Exact SOLUtCLON 
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Periodic Variation in Fin Effectiveness 
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CHAPTER 


INTRODUCTION 


Free convection phenomena occur in a non-isothermal 
fluid under the influence of body force. For example when 
a vertical metal fin is heated in an originally quiescent 
fiuid, heat is conducted! into the’ surrounding fluid to 
cause temperature difference, which in turn generates the 
density variations in the fluid; a boundary layer flow 
adjacent to the fin is then formed as a result of the 
density variations and the body force field, and conse- 
quently more heat is transferred from the fin to the fluid. 

In practical application, stationary fins have been 
frequently employed to increase the total heat transfer 
rate between a solid body and its surrounding fluid of 
different temperature by sincreasing the vareator the contact 
surface. The free convection phenomena associated with 
these stationary fins are normally een to. the runafiiorm 
Gravieational forceytield. 

In this thesis, the heat transfer problem of a 
heated radial: fin rotating simultaneously withivambient: air 
abouteanhorizontal exrsmisastvdied. » Pie pimps clocatedtat 
a certain distance from the center of rotation so that the 


Centrifugal saccelerationuat theviingr oot ican always be set 
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greater than the periodically changing effect of gravita- 
tional force by adjusting the rotating speed. Formulation 
of this problem shows that periodic terms correlated to the 
gravitational acceleration appear in the derived governing 
equations and the solutions may have finite fluctuating 
amplitudes. 

During recent”years” considérabie” interest! has! been 
shown in periodic phenomena associated with fluid flow and 
convection heat transfer [1-25]*. In most of these studies, 
the periodicity with respect to time normally originates 
from a given oscillating boundary condition [1-15]; for 
example, oscillation of a solid wall or its temperature, 
Fluctuation of the free stream fluid velocity or temperature. 
Very few works [16-25] have dealt with problems in which the 
periodicity appears solely in the governing equations. In 
general, itlas Very Gitfri cule to solve the complete set of 
poOverning equations Or these periodic propiems cdrreccty. 
Because of this, several particular Rei bee toe have been 
deve loped. 

The method of successive approximation consists of 
splitting the solutions into two ro more time-dependent 
components, each of which satisfies a different set of equa- 
tions simplified from the original governing equations after 


careful, physical reasoning. Usually, these simplified 


* Numbers designate references given on page 143 
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Scuacrons varesmuch@ecasier Co Solve;s for instance, they may 
consist of synchronous solutions which can be obtained by 
means of separation of variables or Laplace transformation. 
However, this method demands artistically intuitive skill 
on separating the equations and their solutions. Such 
technique becomes very complicated and tedious if higher- 
order approximations are required to produce more accurate 
results. 

The best-known method for solving unsteady problems 
with small fluctuations is the perturbation expansion in 
powers of the oscillation amplitude which has first been 
normalized with respect to the mean amplitude. The zeroth- 
order equations of the expansion are related to the steady 
Dasicet low wiile che higher-order perturbation equations 
are associated with the superimposed fluctuations. The 
accuracy of the solutions obtained can easily be estimated 
PEOM. the Bae Oe perturbation been used. In other words, 
this method is powerful only when the oscillation amp litude 
is small; otherwise work may increase considerably in order 
to include higher-order perturbations for more accurate 
Solutions. | Lt is also very dukely that the method will fairl 
completely if the normalized oscillation amplitude is not 
mich less than unity. 

A powerful technique to handle problems with finite 


c 


amplitude oscillations has been introduced by C.C. Lin [26] 
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who divides the governing equations, as well as their solu- 
tions, into time-average and time-dependent parts. The 
complete set of time-average equations describing the mean 
flow is retained; the time-dependent equations, usually 
linear and uncoupled from each other, are simplified to 
such convenient forms that their eoteeions can be obtained 
analytically. These approximate, time-dependent results 
are then substituted into time-average equations to get 
mean-flow solutions. Hence, the accuracy of the overall 
solutions depends strongly upon the proper simplification 
of the time-dependent equations. This simplification 
generally requires a condition of very high oscillating 
frequency if the oscillating amplitude is relatively large 
in comparison with the mean value. This requirement becomes 
a severe limitation on Lin's method. It appears that no 
rene een using Lin's method can be successfully applied 
to periodic problems under low-frequency, finite-amplitude 
oscillations [ll, 12, 15]. 

In this study, the periodic heat transfer problem 
are governed by a set of time-dependent, nonlinear partial 
differential equations, and none of the three methods above 
is considered to be proper for seeking the solutions under 
low frequency, moderate-amplitude fluctuations. 

Benefiting from the availability of high-speed elec- 


tronic computers numerical methods provide a good possibility 
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for solving such a periodic problem. A consequent challenge 
me the numerical instability of the governing equations 
heightened by the need for unknown initial solutions at an 
undetermined time of the period; a poorly chosen initial 
time level, or an inaccurately guessed initial result can 
easily lead to solution divergence or numerical overflow. 
A proper approach is developed in this study to alleviate 
such difficulties in starting the numerical integration. 

At each time level, the governing equations express 
a distinct, coupled conduction-convection heat transfer 
problem, ana Solution matching fauatne isolvd-fluid anterface 
Mepuececcary Iie ?traislormation developed in a similar 
work [29] is applied so that the coupled governing equations 
can be solved numerically with a "local similarity" technique 
[27-30]. As a result, not only the matching problem of 
interfacial temperature and heat flux at each time level is 
resolved spatially in a few iterative loops, but also the 
periodic conditions of temperature, velocity, and their 
derivatives with respect to either time or any base vector 
are successfully satisfied at every time step through a 


couple of timewise iterations. 
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CHAPTER aie 


THEORETICAL FORMULATION 


The system under consideration is iilustrated 
schematically in Figure 2.1 which indicates a triangular 
fin in a surrounding fluid rotating at constant speed about 
eunoOrl zontal) axis.) Heat 2s) supplvedmat the root of the, fin 
and conducted into the fluid from the fin surface. Asa 
result, a free convection problem induced by a periodic 
gravitational force and a spatially varying centrifugal 
force takes place. 


2.1 BOUNDARY LAYER PROBLEM 


It is assumed that: 

(1), @he £luid as Newtonian. 

(2)- The flow is ‘laminar. 

(3) The flow is two dimensional in space. 

(4) Bela apasiore rts are constant, with the 
exception of fluid density to which the 
Oberback-Boussinesq approximation can be 
applied. 

(5) There is no heat generation in the fluid. 

(6) Wall effect and the trailing-edge effect 
are insignificant. 


(/) Radiation ws insignisican:. 
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CENTER OF ROTATION 
OF SYSTEM 


System of Radial Rotating Fin 


(8) The mechanical work against the body force field 
PS InSLginet Cait. 

(9) The fin slenderness is large and the effect of 

inclination and curvature become negligible [(44} 

(10) The rotation of the system is so fast that the 

centrifugal acceleration is always greater than 
the gravitational acceleration. 

The governing equations of this periodic free- 
convection problem have been derived in Appendix A. Referring 
to equations (A.14)-(A.17) in Appendix A, one can easily 
obtain the following unsteady boundary layer equations (see 


also Appendix B): 


Oa eer (221) 
OS ee eae ee 1 eee a VU - 2ro0V 
t Xx i p.—X YY 
C2582) 
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Sele) met ug COS (Ot) l(b —iieg), 
whee dae Fi 
Te PAU SEN i eS 2d wU 
C73) 
2 
aor ue Cuca TCG) a iced ee), 
and 
Deeley Wie temoaks (2m 4) 
where A=l1 or -1l for the leading and trailing faces, 
respectively. The imposed boundary conditions are: 
X= 0: UD Se Oita ees 
¥ = 0(12xX>0): U=V=0 Tae wee ts) (2505) 
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where T. (x,t) is the interfacial temperature distribution 
of the -Soliderinwand the fadjacentsciuica. 
2.2 CONDUCTION PROBLEM 

The basic notation and the coordinates of the 
PEenduceron neat transter problemeaneche finare ildustrated 
in figure (2.2). It has been shown in previous work [27] 
that) if the*Biot number (Bi) as much less than order one, 
this conduction problem becomes spatially quasi-one-dimensional, 
i.e. the lateral temperature gradient across the fin is 
comparatively small ~A tacit implication. of this small, 
lateral temperature variation is that the thermal boundary 
layers on both sides of the fin are symmetric about the 
X-axis. It is known that this symmetry holds when the 
whole system is spinning at very high angular velocity 
[27,28], and the convective heat transfer coefficients hj (X,t) 
enc hott), shown tn tigures (2.2)) are adenticalseand equal EO 


-K,T, (%,0,t) 


ho(k7t) = 7ho(%,e)) = IDA OCS RS 
Orniwe co 


ASS) 


Moreover, if there is no heat generation within the fin and 
tne Properties Of | cheptin are yconstant, the temperature 
TY (Xt) of this high-speed rotating fin satisfies the 


following unsteady, quasi-one-dimensional conduction equation: 
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BPagure 2.2 Dimensions of Triangular Fin 
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InP this, stuay,, the roOtatzonmspeeda or the system 
is not very high and+h,(X,t) May not. equal h.,(xX,t) unless 
thermal boundary-layer on both sides of the fin are 
symmetric. Further study of the boundary layer equations 
(2.1)-(2.4) is required before the validity of the quasi- 
one-dimensional equation (2.7) can be confirmed. 
2.3 SEPARATION OF EQUATIONS 

A lack of initial conditions is a unique challenge 
in the present problem. The governing equations (2.1)-(2.4) 
are a set of strongly coupled, nonlinear partial differential 
equations with periodic coefficients. The physical conditions 
imposed on the solutions indicate that this is a two dimen- 
Sional, boundary value problem varying periodically with 
time. In principle, a periodic problem like this does 
not necessitate any initial conditions; this substantial 
dviterence from theyordinary, transient-type unsteady problem 
can lead to great GEEPUCUILE yell OobaIning asSOLUtlOn. 
Only in some special cases when the periodic solutions 
are synchronous, i.e. in the product forms of a time- 
dependent and a time-independent functions, does such an 
absence Of initial conditions not*cause any technical 
problem. Unfortunately, synchronous solutions are not 
obtainable in the present problem, and the powerful method 


of similarity transformations can neither remove the time 
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ie 


dependency from the governing equations as well as their 
boundary conditions, nor Simplify the problem by mapping 
it onto a two-dimensional domain. It is thus imperative 
to overcome the difficulty in obtaining proper initial 
conditions so that a successful method ensuring accurate 
solutions can be applied. 

A perturbation method can provide the first degree 
euproximacion’ to the solutions. “When the centrifugal 
acceleration is very large [27, 28], the problem reaches 
quasi-steady state and does not need any initial condition, 
it is logical to start with problems under the influence 
of relatively large centrifugal force field. The solutions . 
are first expanded into a series of linearly independent, 
synchronous functions. Second, the quasi-steady-state 
solutions can be used as the first terms of these series 
expansions. Although in principle, the accuracy of each 
approximate solution can be improved by including more 
synchronous terms in the series, this method is regarded 
as neither reliable nor economical if the rotating speed 
is reduced considerably. Direct integration of the 
governing equations (2.1)-(2.4) by ees of a numerical 
technique is certainly a rational alternative for handling 


such moderate rotating-speed problem. In order to directly 


Py l 


iG 


integrate sequations (271) —(2.4)" vali the missing initial 

conditions have to be guessed at every point on the X-Y 

plane. These guessed values should be accurate enough 

so that after reasonable timewise iterations, not only 

all the dependent variables but also their derivatives 

have to be matched at each point throughout the two 

dimensional domain. Otherwise, all the effort may only 

lead to long computation, solution oscillation and, 

very possibly, numerical overflow, Apparently, the high 

rotating speed solutions, if available, may be a 

Suitable approximation to these guessed initial conditions. 
Although product-form solutions are not obtainable 

in this problem, one can easily split the solutions into 

Steady and periodic parts denoted respectively by super- 


Secepirce seendep sit eO-—l ie substitutions Of 


6 (X,¥,t) 


= 9° (x,y) + oF (x,Y¥,t) 
U(X,Y,t) = U> (X,Y) + UP (x,yY,t) 
5 p (2285) 
V(X,¥,t) = V>)x,¥) + VP OY, 6) 
and eben =) is (Onna) oe SIO Manta 


into equations (2.1)-(2.4), and separations of these 
equations into time-dependent and time-independent parts 


lead’ to 


Ss eS) Soe p 
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= vue + 2dwvP - [Bgcos(wt)]O® - Bw? (R-x) OP 
= [8gcos (at) 16P| (2-210) 
U=vy> + vov> + = p= - wy? = 2iau> + Bu-yo® 
x wy. Poo Y yy (21s 
ea Pa Pele pects DS pS Pee POEs 
Ive SU ec. 1G pe ION be NES i WaNEt Dany ie + DY 
= yee - 2dwuP — [rABgsin(wt)]0° + Bw*yeP - [ABgsin (ot) 10?| 
and 
S,S SoS _ Soe lese Sop Sap Pos Pos 
U oO. + V ee Ke vy lee a oe + V os + U a + V oe 
PoP PoP @P 
+ uPeP + yPor Ke ah (2) 


The right-hand-sides of these equations are functions 
of time and, individually, become zero in the steady case. 
In other words, the steady solutions oae ee Ge and P® 


satisfy the equations 


iS Soke 
Us <2 Vy = 0 CORRES) 
Ss ofevey Ee rll ae: Sen S Boma Ss 
Oks oe VS amex OS 2rWV~ + Bw? (R-X)09 (2A) 
Shes: See oy GS dk ets Ss Sie ens 
Uv, + V Ma pay + Wyy + 2 WU STC mage) (25S) 
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and the boundary conditions 


Kq= tO 0 eo | oo se 0 
Y= Oi (Lex) Ua = 0 6k =e (x) (2ea07)) 
Veo aoe eee Uge = 00 1 Ores ps ag 


This satisfication is always found because of the 
way the steady solutions are defined. Consequently, the 


periodic part of the boundary layer problem is governed by 


equations 


ig ae 
Uy. ar Vy = 0 (2S) 
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+ [Bgcos(wt)]0”, (2.79) 
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The associated boundary conditions are 


Kea Oils RUS = 80 eLOe t=1.0 
race ON <0) ee Ur i= OP en ome tat). (2522) 
Vee aU c= 0) 6 Olt=). pono 


SOolutrvons, Of equations. (2. 03)—(2.16) and boundary 
conditions (2.17) have been studied in previous work [27,28] 
where the conduction problem in the fin is found to be also 
quasi-one-dimensional. As mentioned in section 2.2, one 
of the necessary conditions for a quasi-one-dimensional 
fin equations is the symmetry of the two thermal boundary 
Wayers adjacent to this solid fin. In the study, there’is 
a need for further investigation of the separate boundary 
Payerrequations (2.13) —(2.16) ,and (2.18)—(2. 21)% 
lita NORMALIZED EQUATIONS 


Upon using the normalized variables [27] 


x = X/L y=yY Ral/*/L 
u> = Unie Ral’? vo = Vay ice EW Woe, 
a ae 7a Pennine fe eye 


it has been shown in Appendix B that for air, the quasi- 
steady boundary layer equations (2.13)-(2.16) can be 
simplified through the processes of normalization and 


Orderings | oByY rererringsto equations, (B14) —¢B. 7) in 
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section B.1, these simplified boundary layer equations in 


the normalized form are 


S 
ane vy = 0, (2a24) 
ih ss Ses a ees Ss 
Be (u uy ata Wi u,) = Yoy rr LURE eG iO eae C22) 
Sere Gree 6s 
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and 
s.s s,s s 
a = Ziel), 
u ve Vv vy ae. ( ) 


Bquation (2.26) indicates that the Coriolis force 
merely induces the lateral pressure gradient in this 
quasi-steady part of the problem. Having suppressed the 
pressure term from equation (2.25), one can find the 
pressure distribution p° without any difficulty once the 


: Ss 
SOLucCLON. uu 


is at hand. Hence equation (2.26) can be 
ignored in this quasi-steady problem unless one is interested 
in the pressure profile. Hereinafter, the steady problem 


becomes that of solving the other three partial differential 


equations with the associated boundary conditions: 


Ss Ss 


2 Os uu = 0 go = 0 
V7 Oe Clex 20) ue = vy" = 0 go? = o* (x) (OL, PB) ) 
y=: u® = 0 o° = 0 


Similarly, the periodic boundary layer equations 
(2.18)-(2.21) can be simplified after the following normalized 


variables, obtained in section B.2 , have been introduced: 
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x = X/L y = yRai/* sz T= wt/2T 

uP = UREEr /peRe vP = WELEr/iecRay * oP = (0? /9 _)Fr 
and 

as COL LORE if s<o[ll, 
Ors 

pP t= 2mPPrr/p,,« ,wRal/* iio AOL, 


where parameter 


ah iy, 


( yifs2 
20 Pree 
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Se 


is the magnitude ratio of the Coriolis force tern, 2wu?, to 
the buoyance force term, BgO°sin (wt), PNeeCOMatt On wt. 22 0, )ic 
Physically, this parameter S indicates whether the 2wuP 
or Bg” sin (wt) term gives rise to the pressure /fluctuation 
pP of this periodic problem. 

With reference to equations (B.20b), (B.21b), and 
(B.23b) in Appendix B, the simplified forms of equations 


(2218) =(27 20) are 

= 0, (25299) 
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1_4(,P,P Dey etn P 
+ [ee] (ut oy tv Sl bdyy, (2-32) 
subject to the associated boundary conditions: 
x= a0F i: uP = 0 go” = 0 
Va Orla x0) a: uP = vP = 9 oP = oF (x,T) (2ee32) 
CES uP = 0 t= 0 


As described in section ie., the reciprocal of: the 


parameter A, i.e. 


vee i (Pr Ms 
Ae 25 BO. 


indicates how intensively transient behavior is induced by 
the periodic effect of the gravitational force. 

The lateral equation of motion is, again, "de-coupled" 
from the other boundary layer equations (2.29)-(2.31), and 


the pressure fluctuation can be calculated either from 
Py Sapo yigls ee [=]Ao® Sabin Pee) | (2.53380 
bE GO Lis ox (irom 


Py = [S]2au® + do° sin (27T) (Di esea)) 
Beto Old | er enOtner words; culceper1Oo0rCc portion OL ele 
problem is governed by equations (2.29)-(2.31) and the 
boundary ‘condition (2.52), 118 ONewas NOL. interes codsin =the 
pressure distribution at all. ss 

Equations (2.29)-(2.31) show that the periodic 


thermal boundary layers on either side of the fin are 
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Symmetric to each other. This is consistent with the 
Possibility sthatythe periodic conduction in the Solid fin 
may be quasi-one-dimensional in space. 

Reunion of the split, steady and periodic solutions 


generates the normalized solutions of the overall problen: 


WX, V;T) = u” (x,y) + [ec] uP (x,y,T) 
VAR PY, 1.) i= v> (x,y) + [ool v? (x,y,T) (25934) 
b(x,y,t) = oF (x,y) + [eel oP (x,y,7) 


These solutions satisfy the combined governing equations 
of the boundary layer problem: 


Ls toe ean) (203 5)) 


arll=lu, or uu, <2 vu,,} = Uy + (l-yx)o + [ec] ocos (20T)), 
(275726) 


oe ; 
[glé, tub, + vd =o, Oa 


yy 


and the boundary conditions: 


x = 0: ea 0 d= 0 


y = 0(0<x<1): Us pVetoanU P= G(T) (2.38) 
Vy = @% bl = (0) d= 0 


Again, equations (ees) Oren) indicate the symmetry 
of the two thermal boundary layers; hence the quasi-one- 
dimensional, unsteady fin equation (2.7) may be employed 


Providing; sof COUrSe; the DEOL number a 51) sor athe probplen 
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is small. By using the dimensionless variables 
x = X/L y= Y Rae ye 
and od = CLT) eles leslie 


equation (2.7) can be written as 


36 a7 6 o¢ | 
[el xs—? ani oke ro + aia + [C}$,, (x, 0) = 0 


where 
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NUMERICAt, TECHNIQUES 


PN SeQuatvonsis( 2.36) and iene) peat they transient 


term are expressed by ee ea p 
u, = [spl uy 
(3511) 
tle Pp 
aie Oe etm ee 
where ub and Om are normalized quantities, the boundary 


Wavyerm (equations (2.35) — (22.37) canbe written into 


Ue ey Gare) 


jes) WE or a pied + (l-yx)9¢ 


(3:53) 
vi — pcos (27T), 


and 


a 


Graal Meath ig ae hb (3.4) 


AB 
The first approach is to solve these equations for 
large Fr values. It can be seen from the above equations 
that the fluctuation is introduced into the problem 
through the gravitational force term, $cos(2mt), the 
magnitude of which is controlled by the parameter 1/Fr. 


When the angular velocity w of the system is extremely high, 


Fr = 


implies that the problem becomes quasi-steady [27,28]. 
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in general, the Fr°value is finite and the problem 
becomes periodic and much more difficult to solve. 
Hence, it is natural to seek solutions for very large 
values of Fr in the first place and then progress to 
cases when Fr is of order one. 
3.1 QUASI-STEADY-STATE SOLUTIONS 

It is imperative to develop an efficient method 
of obtaining the high Ss Golimeicins. As described in 
reference [27], this is a Ree ede eens free 
convection problem coupled mien eanottes quasi-one-dimensional 
conduction problem through their interfacial boundary 
conditions; the numerical solutions are obtainable through 
the iterative matchings of both temperature and heat flux 
at the solid-fluid interface. Although use of the local 
Similarity transformation and the numerical technique 
employed in previous work (27,281 has assured this solution 
matching, an improved scheme to provide faster numerical 
integration as well as solution convergence is still 
necessary because in solving the low rotating speed problems, 
numerical computation may increase enormously in seeking 
the periodic solution through timewise iterations. The 
development of a faster scheme in obtaining solution at 
each time step can drastically reduce the computation time. 


To this end, a slightly different method is described below. 
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3.1.1 BOUNDARY LAYER EQUATIONS 
It is convenient to transform the governing equations 
keine 4) 2 125') Pande (2027) hinteganothemiset of partial 


differential equations by putting [29] 


a epee Lael, Sejace aly 
co ne (PEARY 
(3 iD) 
Ss Mowe L evry Ge | 4E\ sy wes 
Wo (x,y) = ces) (3): Se eo) 
and S Ss 
) (x,y) = © (oan) 
where v> (x,y) is a stream function, defined by 
ue = vy and ; vo = ve (3%.6)) 


such that equation (2.24) is automatically satisfied. The 


newly transformed equations are 


Ss Ss-S 2 S. Ss 
P + - = + + co) 
5e eee f£ fon 3 (£7) Clr ss) 
| (Siw) 
43, S28 Ss .s RPS 
= = £ ren + (1+Pr) ¢ 
lethal eer) O) 
cS} s,s 4 s.s Sia 
+ = - Oo £f Cre: 
subject to the associated boundary conditions. 
TpseOe W<EAWE ba emiche ci Matendal =) 52) 
G3 5.9)) 
Sea eh f° = 0, or =O. 


It is worthwhile mentioning that transformation 
(3.5) is a combination of the Lefevre type [31] and the 
Blasius-Howarth type [32-34] transformations. The Blasius- 


Howarth transformation stretches the coordinates [45] in such 


jing om © 
104 o9 
7 44 ' 
1 wip@ 
od 
Leave 
ral 
at ce 
S57 nee 
a D 
: ae bck Au 7 Sa 
y, ; 


7 


ut ey mt NA 
in : aa ney oe 


. 
i ee 


M mh 7 and . he True ) 


{ i jer ) sarod idsi hecate nip ae WeOrsey oe azeKe 
ten : ey pee - 7 . \ 


rr 


Lt 


— yas) 


| hye 


Se oe, ey uke bey i iteabeee 


é ; F eo ; ah 7 a 
sige! Ana etaeibenlauitsatas 


; = * ; ek: bs, 
¥ « + ‘ny Atal he ‘al a) 
bi Ga4h\-Sdete Py eee awe 
bem : Me Se CaaS? 


a eee es) Ane By hte, A * _ 


Biles a Bric fact satpro Lod ik et 
ro ~ ~*~ * dl ¥ yeu : 
(ay Be: Ba.” fl ae le He stab ¢ 
Pil. =e a : yale i ee uy i NB i as 
aL ‘ . i} F : ne ie ave pai ae 
Berta dine ei ab i, od Lotta Be ae 
ia by iva ‘i ee n : ‘ I - ; : ae 


fs re 


(1 Aes aise apbte al a ve tedtannes ae ae 
iF jac 


: ae ey Cl fas a rea ‘al 


r ‘eh oe i Wy ets 


a way that the ¢-directional dependency of the problem is 
Gieatly reduced, and) hencevalltheys-derivatives of the 
solutions can be replaced more accurately by finite 
difference approximations. The transformation of Lefevre 
SerVveseuie pPuEpose (Of reducing ithe eL rece, of (the Prandtl 
number on the obtained results [35] and thus Pere 

the need to determine separate solutions for a large 
range of fluids. 

The integrating evens of this problem is very 
Eamilarpto that of the: previous work | 27-29). (First, *the 
Wiole;region is discretized, along theve-divection) into 
Paratee lace Levels.) | /SeCOUud, tie «dioCreteyva les Vor che 
unknown interfacial temperature o- (5) are guessed at each 
E-level. Physically, this set of guessed values must 
Picrease MONOLOnLCa lly rromr tne: tip, of, the’ fins Third, “by 
use of a fourth-order Runge-Kutta method equations (3.7) 
and (3.8) are integrated in n-direction at the first €-level, 
where the right hand sides of these equations have been 
suppressed because of €=0, Fourth, the integration at each 
subsequent €-level is carried out laterally across the 
boundary after the §-directional derivatives in equations 
(3.7) and (3.8) have been replaced by finite difference 
approximations. Finally, the discrete values of interfacial 
heat flux, corresponding to the guessed interfacial tempera- 


ture O08), are obtained. 
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An initial-value technique was used in combination 
with the Runge-Kutta method to integrate equations (3.7) 
euavel NGSTan ely oie) Paseo tone Experience showed that numerical 
integration of these boundary layer equations was very 
sensitive to the missing initial values of cs and a at n=0; 
an inaccurate initial guess could easily lead to solution 
overflow. This numerical instability was overcome by the 
continuation method [36] with a least-squares convergence 
criterion on satisfying the asymptotic boundary conditions 
at the edge of the boundary layer (371 - 

Numerical ex ,eriment was conducted to determine 
the step size An for the above Mevadaclnen. By varying 
An from 0.5 to 0.005, integrated solutions were obtained 
and compared; it is found that as An reduced to 0.005, both 
fant 30) and 0° (E,0) converged asymptotically to two 
corresponding values. A step ohipes We An=0.04 was then used 
in the integration so that the accuracies of Ee ite 70) and 
O° (E10) were within one percent of their asymptotic values. 
So ie2 ) FIN, EQUATION 

Equations (2.24)-(2.27) shows that the steady 
thermal boundary layers on both sides of the fin are 
symmetric with each other; hence, the conduction problem 
becomes quasi-one-dimensional when the Biot number (Bi) 


is much less than one. 
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Introducing the normalized variables of the 


boundary layer problem: 


x= X/G, y = yRa!/*/z, t= wt/2r 


and Oye S ea AGI OY) (Us 


into equation (2.7) and separating the resulting dimension- 


less equation after the substitution of 


b(x,y, 0 = 8 (x,y) + [z]gP y,t), 


and 2 1 
Gy batx)e pel og (rt) 


one obtains the quasi-steady fin equation: 


Pal ofan Nc : 
< eo FP APL ae [C] ¢, (x, 0) = 0 eee 0)) 


and the periodic fin equation: 


ie aoe 
[S]x5—— + xX ae ae f + Lele, G7) = Q, GSh ets) 


ax4 


The parameter C, already defined in previous work [27] as 


K 
pa ee pei 
Somer areas 


’ 
is a determining factor of the coupling between conduction 
and convection systems; the new parameter B defined by 

B = 20K. /wL* 


indicates the insensitivity of the fin temperature in 
response sco the ineat flux, Cluctuation imposedt ony tie, fin 


surface. 
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Bya vl rec .OLmtheseranst ormacroumus> jy, equation 


(a LODs Canwoe Written: as 


a265 ao® 
E : + sE- fi Regen ae) = 0 vice 
ag 
Ons Ss 
a 7? Be 
Beare) + RSE, OF = 0 (3.13) 


subject to the boundary conditions 


Ce sary = 25 = 1 
S 
Saco OR pita = Shes (3.14) 


Sunes ane fin SA, Miner in He. form of (3.12) ae 
(Srl so) eroustigularcvatec=O;.010 2) CULEVnOCCULS JGuring): Lts 
numerical integration. To accommodate this, one may use 
the following approach. 
If the functions F(E) and H(&) are defined, 
respectively, as 
gta sass 
R(E)) = qe \Sapre ands, | H(e ji ean ’ 
equation (3.13) can be integrated from €=y to €=A& with 
a fourth-order Runge-Kutta method after the missing boundary 
condition H(yY) is guessed. Once the discrete values of F(&) 


and H(é) at each § level other than §=0 are at hand, one 


Gan icalculate H(0) with theaid of the: extropolation formula 
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[46] : 

H(O) = H(A€) - So [ S5F (AE) - 59F(2AE) 

+> 37E (SNE) = OF (ARE) IPs OL KE”) , 

and then modify H(y) so that the boundary condition H(0)=0 
is satisfied. The evaluation of F(§) in the above formula 
Uomoimplesscince 7equacion (3.13) indreares that © (6) 1s a 
function of neither H(&) nor Bee Consequently, the 
Mrcegrat OL this ic), 1.6. He}, 1svalso nothing more’ than 
a function of §°and the known function 7 (E,0). The value 
of H(0) calculated from formula above is directly related 
to the initial guess H(yY); the exact value of H(y), which 
leads to the satisfication of H(0)=0, can then be found from 


H(y)-H(0). With AE<E<y, de® HE) 
do iea as 


arevOnVvlOusly4~ODtLatnable; sruLrchery tn cregravion gives a new 
set of discrete temperatures oor among which the value 


of 0° (0) is, againvoocained) through exctropolation.) That 1s, 


dom (Az) doe (2A) 
e5(0) = a (xe) = 4 [Se = ee 
47to (348) pita 408) 
OT nara api! eaated 
=P Ge) 1558 (ag) = 2? (2a) 
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Bes) MATCHING 

The whole problem of this steady heat transfer 
model has been shown to consist of two sub-problems 
coupled together through their interfacial temperature and 
interfacial heat flux. The two-dimensional free-convection 
problem is governed by two boundary layer equations (3.7) 
and (3.8); the quasi-one-dimensional conduction problem 
Homieserlped bythe rim equation, (3-02) 0r 1(3 013)... The 
matching of the interfacial heat flux has been considered 
wiwcne LOrmulatiom-or equations, (3.62) and (3.13) >) hence, 
the only interfacial matching left to be done is on the 
temperature. In other words, the assumed wall temperature 
0° (E) of the boundary layer problem must satisfy the 
conduction equation of the fin and its associated boundary 
conditions. 

An iteration method is applied to obtain the 
interfacial temperature profile OLE One hte hrs i, 6° (= ,0) 
is guessed and treated as a known boundary condition of 
the boundary layer equations; the integration of these 
equations can, in turn, provide the unknown interfacial 
heat flux 7 (&,0) corresponding to the guessed SE, 0)". 

In distinction to previous work [27,28], only the obtained 
cents) is fed into the fin equation to generate a new 


interfacial temperature distribution o°(E). The guessed 
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6° (E,0) is then compared with this 0° (E) and modification 
of the former is required whenever a difference greater 
‘ase lune see detected. With the newly modified 6°(&,0), 
the procedure is repeated in the hope that a new 0° (8) 
can match this be 80") to a better degree. Hence, this 
interfacial matching problem technically becomes a matter 
of finding an efficient way to improve the initial Hise 
Dae, Oye 

Practical experience indicates that if the para- 
meter C of equation (3511) is teas onan one, ecie 
average value of the previous a” (E,0) and thercalcutated 
O° (é) is a good initial guess of oT CE sOD for the next 
iteration loop; generally, the interfacial matching is 
achievable ina few iterations. This is not true, however, 
when the parameter C is greater than One. On the 
contrary, this modification method by averaging may some- 
tiNesmicad LO tile soscrt atrons on even) divergence Of tne 
interfacial temperature. TOyaccommodate this; this averaging 
method can only be used in the first iteration loop and 
after that, a linear or nonlinear optimisation technique 
[38] should be employed in getting the improved initial 
guess one 0) It has been found that in this way the 
convergent criterion of the interfacial temperature profile 
can be satisfied to an accuracy of Ghats within 3-10 iteration 


loops. 
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Thescesult of spreviousyworwe27.25) pecan suprovide 
good initial guess for Dente) One Although the local 
Similarity transformations of the present and previous 
Studies are different, a closer inspection of these 
transformations discloses that the interfacial temperature 
DLOCULew Sa tunCceLonvoT theasame, Giviariabwle, jan et cases 
and, hence, does not vary with the transformation. For 
this reason, it is a great advantage to use the existing 
interfacial temperature of previous work as the initial 
guess in the present study so that the converged solution 
can be obtained almost immediately. Another advantage 
which stems from previous observations is that in the non- 
dimensionalized domain (x,y), the converged interfacial 
temperature profiles are insensitive to the variation of 
parameter y. In other words, once the interfacial tempera- 
Pune wcorresponding tol asmalliy value is obtained, it .can 
be used as a good initial guess for problems of the same 
Parameter C but largeriys Values.s s Both advantages: are, very 
favourable as far as the computing time is concerned. 

Several numerical experiments had been carried out 
before the present computer program was thoroughly applied 
to this study. A linear interfacial temperature was chosen 
as the initially guessed input profile to a testing case 
Ofev—OsyeC—l.0 mand Pr-0. 7/2 0 bacetlentvagreenent between 


the obtained solution and the previous result [27,28] was 
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Oba nedewmrotcpaenzesroL, 0.05 and (0c025ewere» chosen for the 
€ variable when y=0.5 which was the greatest y values used 
Pptiiswotudy ., sCOMpertson Of /Ene results isnowed that not 
Only the individual values of tworsolutrzons but also their 
LOce builecatatransrter (rates didpnotvamttersby more «than, 0.5% 
from each other. To save computer time and storage, eleven 
S-sitecosc COLresponding tovAs=0.05) were used Lor. y=0.5. 


For convenience, same number of &-levels was used 


throughout the study even for the cases of smaller y values. 


S225 (PERLODIC SOLUTIONS 

Two methods can be used to solve this periodic 
problem pending on the fluctuating magnitude of the 
solutions. The boundary layer equation (2.36) indicates 
that the periodic source of this problem comes from the 
gravitational force term, ?Cos(27T), which increases its 
effect with the parameter 1/Fr. If this parameter is much 
less than unity, a perturbation-type series expansion is 
a proper method of simplifying the governing equations 
before the use of the numerical integration technique. 
Tretwis parameter el /Preicanotemucny tess uthan Unrry,) the 
perturbation method becomes tedious and uneconomical 
because of the necessary inclusion of more series terms; 
a direct fapproach of solving the unsteady problem is then 
inevitable. In other words, the periodic solutions are 
divided into two categories: those of small amplitudes 


and those of moderate amplitudes. 
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So2..) OMALL AMPLITUDE 


RELorring) co) transrormativonms(¢se5) pone can intro— 


duce, a local simigarity transformation | 29 | 


ey ee fess) lenis 


Bs Po PA ENS [sh 
(aes) (3) £(E,n) (3.15) 


w(x,y) 


O(x,y) = O(F,n), 
where (x,y) is the stream function defined by 

Day and ve “Vx, (37, 16,) 
into the boundary layer equations (2.35)-(2.37). The 
continuity equation (2.35) is then satisfied completely by 
definition (3.16); the other two boundary layer equations 


(2.36) and (2.37) are respectively converted into 
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Response parameters Ky and K, represent, respectively, the 
Significance of the transient responses in the fluid and 
solid to the periodic effect imposed on the system, while 
parameter K, indicates the degree of coupling between 
conduction and free convection. Transformation of the 


poundary conditions, (2.38) leads: to 


n =-0(0<&<y): fetal 20", 2(€,0,T) = % (6,7) 
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The associated boundary conditions of equation (3.19) are 
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If 1/Fr is small the solutions may be written in the 


perturbation series of the form (Appendix D): 


£(E,n,t) = £°(E,n) + el (E,n)cos (207) 
+ #(E,n)sin(2nt)] + o[—] 
Ere 
6(E,n,t) = 88 (E,n) + S19 (E,n)cos (2mT) 
; ! 1 
+ (E,n)sin(2mt)] + o[——] (B22) 
Er 


where £°(E,n) and 6° (En) are the quasi-steady solutions 
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and % etc. satisfy the simultaneous equations: 
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and the associated boundary conditions: 
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Higher order corrections can be obtained in a similar manner 
(Appendix D) but are not considered here. 

Periodic solutions of large amplitudes (1/Fr>1) 
excluded from this study due to the fact that a different 
physical problem can be expected if the periodic, gravita- 
PLOvaly LOCC mEGrMN (cos Ant) /Preprevalts Over the centrirugal 
ROeCCBCeIM Cl =vx) gem equation) (230) this may berzexplained 
by considering a typical case of Fr=l. At the moment the 
bocacing@tiansds invits) vertically downward. (7=0). postion, 
aforesaid two forces act in the same direction and promote 
a boundary layer flow with its leading edge located at the 
fin tip; however when the fin rotates to its vertically 
upward position (rt=1/2), these two body forces act against 
Geach other and the prevailing gravitational force tends to 
develop a backward flow with a new leading edge located at 
the fin root. This alternating leading edge problem is 
discussed further as a topic for future research in 
section Sh Pe | 

New Grviticuityuinwsolvangvequartions (3. 23)'— (3.28) 
is encountered after the domain (&€,n) is discretized into 
eleven €-levels. Comparison between the quasi-steady and 
perturbation problems indicates that not only the number 
of governing equations but also the number of unknown 


interfacial boundary conditions is doubled in the latter. 
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Tiageel Sy sene witegrationsof boundary, tayer equations (3.23) - 
(3.26) at each €-level requires the surmise of four 
eigenvalues: £, (6,0), £,, (6,0), ®, (6,0), and §(E,0). 

Most of all, the doubling of the interfacial matching 
requirements between two conduction and four convection 
equations is the main source of new difficulty. 

Gwing tothe linearity, of thevperturbation equations 
(3.23)-(3.26), the search for the four eigenvalues at each 
E-level was accomplished within one iteration by using 
the method of adjoints [36]. 

The overall problem was successfully solved by 
means of an iterative scheme. At first, two interfacial 
temperature profiles | and S| which satisfied the 
perturbed boundary conditions of the fin equation were 
assumed. Based on these temperatures, equations (3.23)- 
(3.26) were integrated simultaneously at each —-level to 
yield two temperature gradients o. and Oa at. the solid- 
fluid interface. These temperature gradients were then 
used in the integration of the coupled conductive fin 
equations (3.27) ‘and (3.28) to give two new interfacial 
temperature profiles. If these newly obtained temperatures 
etdawnotmagree with their previously assumed values to LOpe: 
their mean values were then used as refined values. The 
whole procedure was repeated until the convergence 
criterion of THOME for both interfacial temperatures had 
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SB s2i2) MODERATE AMPLITUDE 

A direct numerical approach is used to solve the 
problems of the Froude number being not much greater 
than one. Reconsideration is given to the following 


governing equations: 
pap 
[Deve se a = 3f +' (1+Pr) {1+ cos (271T) /Fr] o 


= ER(f foe - fppfe) + (LPr) 8] + (K,/2m)EY%e) (3.17) 


j | 

On + £9, = gelf,o, - Hf.) + (Ky/2m E29, (3.18) 
3. Sein ioe ty) 

Bea +o KE" (E,07T) = 0 Bicatp, 


eZ eee tl. 


and their boundary conditions: 


96 
ee WE St 
GE = Ys fe = 1 | 
i= 0: £720" f= 0, GPE Ope eee) 
n= © f= 0, #=0 


Now, if one knows the exact values of time-derivatives 
eS and on ate everypDOLnte(e,niandmacearspeckritc rime TA 
these equations pose a two dimensional problem which is 
muchveasier to solve. | Once the” solutvons Jat this jnitiral 


time level have been obtained, the time-derivatives Be and 


ae at the second time level can be approximated by backward 
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finite differences. Again, a two-dimensional problem is 
produced. Similar approximations can be conducted in the 
Successive levels over a whole period of time. 

THe Mal WAL ee VCwLEy Or ethie waporoaan is yin) the 
pUudrcCrLOUuSHChoLce of thewinitial timesderivatives eT 
and oe at every point (&,n). This choice should avoid 
causing any numerical instability during follow-up 
integration, and should lead to fast convergence of the 
solutions. A consequent problem is to determine the 
proper time level for T; SO that higher accuracy in approxi- 
mating those transient terms Sa and a may be achieved 
at least at this initial level, and even better if it can 
also be done at subsequent time levels. 

When the Froude number is moderate and much greater 
than one, previously obtained perturbation solutions are 
relatively reliable and hence provide good approximations 
to the initial values of Te and o- The remaining problem 
is to determine a proper initial time so that good approxi- 
Mations can also be obtained at subsequent time levels. 

At the second time level, the accuracy of the two- 
point finite difference formula in approximating those 
time derivatives Aeoetie strongly upon the product of the 
time interval At and the magnitude of the second order 


time-derivatives. Because of the long computation required 
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in solving the two-dimensional problem at each tlevel, the 
time interval At should be as large as possible. Hence, 

the onlyjaiternative to obtain*higher accuracy in using 

the two-point finite difference approximation is to start 
solving the problem at the specific time level, after 

which the perturbation method detects the minimum magnitudes 
in the second order time-derivatives. 

Owing to the fact that the total heat transfer 
rate reflects the overall features of the problem, the 
time level where the total heat transfer rate of pertur- 
Dacion result sndicates a=straignt); line, ora) point of 
inflexion with respect to time, is an appropriate point 
for the second time level if Fr is slightly greater than 
order one. In other words, the appropriate initial time 
level is just the level one time interval ahead of this 
point. 

When Fr is not much greater than unity, the pertur- 
bation results are unreliable; indiscriminate use of them 
in either approximating initial transient terms or deter- 
mining the initial time level can easily lead to numerical 
overflow. However, a proper decision may be deduced from 
careful comparison of various manifestations among: 
(Pyethesperturbatvomeresulu rl otminewlarger hr ocase, 

(2) the predicted solution sl of Pl, where Sl is obtained 
through the direct integration of periodic equations 
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(3)) the perturbation result (P2of the present “(small Fr) 
problem. 

The total heat transfer rate Q can be used in 
selecting the initial’ time level. Lt is’ logical \to assume 
tiatevartation) in the value,of Olissinaphase with the 
overall changing trends of system temperature and velocity. 
A zero variation in Q indicates a momentarily unchanged 
system which, in turn, corresponds to zero or extremely 
small transient terms in equations (3.17)-(3.19). 

If comparison among solutions Pl, P2, and Sl does not 
reveal severe or irregular phase shifts, it is reasonable 
to choose initial time level from those of maximum or 
Minimum Q values so that the inevitable errors in approxi- 
Mating transient terms become less effective because of 
their comparatively small magnitudes in the equations. 
Several most typical cases are illustrated in the following 
paragraphs and the guideline on determination of initial 
time level will be reinforced by these examples as the 
discussion proceeds. 

Figure (3.1) shows one of the most favorable cases. 
The perturbation result Pl for a large Froude number is in 
phase with its Preaiened SOLUGLON Suen addition, the 
perturbation result P2 for a moderate Froude number does 
not reveal significant phase-shift from either Pl or Sl. 

Tn this case, the solutionyS2ypredverecdupype2 for the 


moderate Froude number can be expected, with more confidence, 
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COsbe inepnase with all’ the -solutgionsman Figure, (3.1). 
More specifically, both leyels of t=1/4 and t=3/4 can be 
chosen as t;. However, the t=1/4 level is considered to 
be proper because of its smaller magnitude-difference 
between Pl and Sl. Once this initial time level has been 
determined, the transient terms in equations (3.17)-(3.19) 
can be approximated by using its perturbation solution P2. 
On the other hand, one can choose T=0 or 1/2 as Ts and use 
SimerOrmtne i naetLalssolutions )2t higher accuracy in the 
initial finite difference approximation is more desirable. 
Sometimes, the phase of the solution may be 
sensitive to the variation of Froude number, as is shown 
in Figure (3.2). Although the perturbation result Pl 
for a large Froude number is in phase with its corresponding 
Soluttonesl) ps Ceasecsignitrcantiy OUL OL Dhaseyrrom.P2. 
Since there is no phase=shirt between SI) and Pl) it is 
plausible to expect that no phase-shift between P2 and 
S2, and both t=0 and t=1/2 can be the candidates for Tse 
However, t=0 is properer because of its smaller curvature. 
The reason behind this is: that when the curvature of the 
selected extreme point is small, the error in approximating 
those transient term of governing equations (3.17)-(3.19) 


with P2 becomes less susceptible to any unpredeicted, but 


possible, phase-shift between P2 and S2. In addition, the 
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second order time-derivatives are proportional to the 
curvature and a smaller curvature means the possibility 
of smaller truncation errors in the use of two-point 
finite-difference approximation at the second time-level. 

Concerning the example of Figure (3.3), the 
perturbation approximation Pl is about one-eighth of a 
period out of phase from its predicted solution Sl. 
Fortunately, this Pl curve is in phase with the pertur- 
MaLlonasO MULLOn) P2. If) Pl and P2paremnot very Lar apart 
in their magnitudes of Froude number, one can assume that 
the Froude number does not have a significant phase-shift 
effect to the solutions either. In contrast with two 
previous examples, Sl in lieu of P2 is more suitable for 
determining the initial time level as well as for 
approximating all the transient terms of the governing 
equations (3.17)-(3.19). Again, the level of t=3/8 rather 
than’ t=7/8 is a better -<choice for T; because of its 
smaller curvature. 

A more difficult case is illustrated in Figure 
(3.4), where the perturbation solution Pl is in phase with 
neither its precicted result Sl nor the perturbation solution 
P2 corresponding to a smaller Fr. However, comaprison 
between Sl and Pl shows that the extreme points of the 


former locate about one-eithth of the period behind their 
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corresponding points of the latter. “Assuming the same 
correlation also exists between P2 and its predicted 
solution S2 for a small Froude number, one can choose 
the level of t=5/8 as na and start solving equations 
oe l7=\(3019)) by approximating all the transient terms 
with P2. 

Sometimes, the change of Froude atbayrc can cause 
a different phase-shift in the solutions. The discussion 
of the four examples above is based on the Bee Con that 
the correlation observed between solution Sl and its 
approximation Pl corresponding to a large Froude number 
also holds for smaller Froude numbers. This is not always 
true and the correlation may vary with the Froude number. 
As depicted in Figure (3.5), the solution Sl corresponding 
to a large Froude number has a time-lag of Ll behind its 
perturbation approximation Pl, whereas the unknown solution 
S2 for a small Froude number may have a double time-lag 
Or L2 behind its perturbation approximation: P2. | Apparently 
this is a good example of poor accuracy in perturbation 
solution when the Froude number is not large. If one 
wrongly chooses tT; ~level to be 1r=3/8 instead of 1t=1/2, the 
numerical integration may eventually lead to divergence of 
the solution. To overcome this, a proper correlation 
between the solution phase-shift and the variation of 


Froude number should be established by observing additional 
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Solutions corresponding to one or two intermediate Froude 
numbers. 

Figure (3.6) indicates a severe case where not only 
the perturbation result Pl corresponding to a large Froude 
number deviates significantly in magnitude from its approxi- 
mMated solution Sl, but aliso no. correlation exists between 
e-remespoints Of Plvand SM. Moreover. comparing Pl with 
the perturbation solution P2 corresponding to a smaller 
Froude number, one notices that no consistent phase-shift 
holds between their extreme points. That is, the maximum 
point of P2 is one-eighth of a period behind its corresponding 
point (at t=1/4) of Pl while the minimum point of P2 is one 
quarter of a period behind its corresponding point (at 1T=3/4) 
of Pl. It is recommended that a series of problem with 
gradually decreasing Froude number be solved. Owing to the 
lack of accuracy in the perturbation approximation, one 
should use Sl instead of P2 in determining the initial 
time level, as well as in approximating those transient 
terms of equations (3.17)-(3.19) corresponding to the 
new but slightly smaller Froude number. In addition, 
thestime evel of t=-3/4)in) Figure 36) sis) an ‘appropriate 
CGhoLeestoxr Ts because of its smallest curvature among all 
the extreme points of the Sl curve. 


Once the initial time level Ts and its associated 
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time-derivatives one and one are properly determined, 
Pec Gc de lOmeor Cle governing tequacl Onomdt tits eini tials and 
subsequent time levels can be carried out without any 
difficulty. At the end of one time period, the timewise 
convergence of solutions is checked. Theoretically, 

this convergence can be checked at each time level after 
the first timewise iteration, provided that the solutions 
obtained at each time level have been previously stored. 
However, the immense size of computer memory required 

for storing such not-yet-converged solutions makes oo 
impractical to check the convergence of the solutions at 
every time level. In this work, the check was performed 
only at the end of each period and the timewise convergence 


was achieved within 3 iterations to an accuracy of lowe. 
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RESULTS AND DISCUSSION 


As demonstrated in the previous chapter, the results 
of this periodic heat transfer problem can be divided into 
quasi-steady solutions, perturbation approximations, and 
"exact" solutions obtained from the direct integration of 
the governing equations (3.17)-(3.19). These equations 
and their associated boundary conditions indicate that the 
Disb eme contains, Six, parameters: (Pr Er, Y, Ki Koy and 


Gv (eke K3). Among them, the effect of Prandtl number (Pr) 


on the solutions is removed by use of the Lefevre-type 
transformation (Appendix C); hence a fixed value of 0.72 
is used (as the Prandtl number of air) throughout this 
study, and the number of parameters is therefore reduced 
we) (Saaky eave 

The Froude number Fr is the most important 
parameter in this study. Since solving this coupled 
convection-conduction heat transfer problem with the 
necessity of additional time-wise iterations requires 
expensive numerical work even for one set of parameter 
values, it is uneconomical to study the effect of all the 
remaining five parameters; only the parameter of greatest 
interest in the present study should be considered. One 
of the main objectives of the present work is the develop- 


ment of a technique for handling periodic problems of 
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moderated fluctuating amplitude, which strictly depends 
Onwthemmagneeude ofl Er) Wi For thisvreason,salthough some 
attention is also given to the’study of other parameters, 


the investigation is mostly centered around the Froude 


number. 
4.1 QUASI-STEADY RESULTS 

The quasi-steady results can provide the basic 
information to the problem. Although this quasi-steady 
case has been previously investigated [27,28], the present 
results are obtained in a new coordinate system using 
soem transformations. In addition, these quasi-steady 
solutions are the base in constructing perturbation approxi- 
mation which, in turn, may be the key to the solutions of 
the periodic equations (3.17)-(3.19). Thus, presentation 
of these new quasi-steady results is not only for the 
purpose of comparison but also. for the sake of completeness. 
For convenience, these quasi-steady solutions are divided 
into three parts. 

In order that the present results can be compared 
with those of previous work, parameter C rather than K, is 
usedshere.. This is) certainly. acceptable because. for a 


fixed, Prandtlanumberyol On 2, parameter K3 expressed by 


bi 3 Pr —3 1/4 
Revie aera ae be 
Tormeinelact ~san DuncLlon aor. paranecte: ogc sandey amin tis 


study, investigation into the effect of these two parameters 
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Pomcal Piet sOuL with six valuesvor tem uasal mney 1.0. 340 
SOgaidelOno) ands five values ofa (emmenusoy 0.3, 0.4 and 
Ci Sys 
4.1.1 SOLUTIONS OF BOUNDARY LAYER PROBLEM 

As far as the study of this free convection 
problem is concerned, the investigation encompasses the 
solutions within the momentum and thermal boundary layers. 
The shapes of these boundary layers are sketched as dashed 
lines in figure (2.1); namely, the boundary layer thickness 
increases from zero at the fin tip (x=0) to finite value 
at the fin root (x=1). However the definition of new 
coordinate variable n of transformation (3.5) indicates that 
in the transformed coordinates (&,n), the closer the fluid 
flow to the fin tip the greater the stretch of boundary 
layer in n-direction. A proper transformation should 
therefore stretch the boundary layer in such a way that 
its thicknesses at different €-levels are of the same 
order of magnitude. These boundary layer thicknesses 
become independent of € only if similarity solutions exist. 
The possible existence of similarity solutions for this 
boundary layer problem is certainly worth considering. 
4.1.1-1 DEPARTURE FROM .SIMILARITY — 

Figure (4.1) shows the velocity and temperature 
profiles at three € levels along a high-speed fin with 


parameters C and y respectively equal to 5.0 and 0.5; 
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Figure 4.1 Departure of Quasi-Steady Boundary Layer 
Solutions, From,cimimvara ty 
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apparemely, departure from’ similarityers quite: significant 
inebotirspcoliles ofthis boundarywilayer Brom lent However, 
the fact that both thermal and momentum boundary layer 
thicknesses remain almost unchanged along the fin direction 
in the €-n domain, substantiates the coordinate stretching 
feature of the Blasius-Howarth type ee erormstion (BRIS). 
4.1.1-2 EFFECT OF PARAMETER y 

To study this effect, velocity and temperature 
PEOtiles Bacythesroot of the rotatingifin with parameter 
C=1.0 are plotted in figure (4.2) for three values of 
vom(0.17) 0-8 and 025). The velocity varies considerably 
with y, but the temperature changes only slightly. Moreover, 
the thicknesses of the thermal and momentum boundary layers 
are nearly unaffected by y. 
Ae le lof oe Pel, Ob PRRAMETER.C 

Since }the magna tude of parameter C represents) the 
degree of coupling between convection and conduction, inves- 
tigation of this parameter effect seneas Boer order-of- 
Magnitucdes or Ca(0n0 1,0. lb. 0 Pandel0)..0) \forvy=—0'.. 1, as 
shown in figure (4.3). Apparently, parameter C does produce 
a significant effect on both velocity and temperature profiles 
at the root of thersfin. In addition, the effect'on the former 
DEOTileswis bY, faremore) noticeabliemtoan ongtie watter in 
contrast to these, neither the thermal nor the momentum 
boundary layer thickness at this fin root varies noticeably 


with the parameter C. 
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Effect of y on Quasi-Steady Boundary 
Layer Solutions 
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C=.01 


Figure 4.3 Effect of C on Quasi-Steady Boundary 
Layer Solutions 
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4.1.2 SOLUTIONS OF CONDUCTION PROBLEM 
The Main interest of this investigation is'on the 
temperature distribution along the’fin. Yin most studies 
Of the vee ‘convection problems [39-43] heat conduction 
in the solid is ignored, and either a power law or an 
exponential temperature distribution is assumed along the 
SsOlid-flurd-interftace. It is logical to’ inspect the 
validity of these assumptions for such a coupled conduction- 
convection problem. 
4.1.2-1 EFFECT OF PARAMETER v6 
Figure (4.4) shows the interfacial temperature 

evong —ctueetin Lor | hoursvalues! or Cy(0c0Us 0.) il. 0 sand 10.0) 
and two values of y (0.1 and 0.5). The actual investiga- 
tron wneludes)five values/ ory; (0.1,"10, 25m"0 23, 10 14" and) 0.5) 
for each value of C. All temperature profiles of y=0.2, 
Ono andl0.4%fall between that of /y=0.0 and’ 0.5. For) the 
sake of clearness, only profiles of y=0.1 and 0.5 are 
presented on this plot. 

7 The insensitivity of normalized fin temperature 
to the effect of y provides a favourable feature from the 
economic point of view; that is, once a matched interfacial 
temperature for a problem with specific values of C and 
y is obtained, it can be used as an extremely good initial 


guess of the interfacial temperature for all other problems 
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Interfacial Temperature 
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Of the Same Cy value but,with different parameter y. Large 
amounts of computation are hence saved due to the fast con- 
vergence of the interfacial temperature. 
4.1.2-2 EFFECT OF PARAMETER C 
Figure (4.4) also shows the significant effect of 
Cron, the solid £in temperature. “It as clear that an 
isothermal surface temperature is a good assumption to 
the boundary condition of a free convection problem only 
if the parameter C is very small. As the value of C 
increases, the assumption of constant fin temperature 
becomes increasingly invalid; but for C=o[1], a linear 
profile with a non-zero temperature at the tip of the fin 
is still the character of the fin temperature. The assump- 
tion of zero tip temperature for the fin is only a valid 
approximation when the magnitude of C is much greater than 
Onder one. 
4,1.2-3 COMPARISON WITH POWER LAW ASSUMPTION 
Peele wlaceNedeim ce. hac hau temperature 5 follows 
the conventional power law assumption ge=ax of free con- 
vection problem, it should have the following features: 
i) 9°=0 at x=0, except for the isothermal case where 
n=0.5 “According to figure (474) Vethis as, only possibie 
whemnic> >O;)%. 
a einem Lot lo. $0 vs. x on a log-log scale should reveal 


aus trarghnt: ane 
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To check the yalidity of this power law assumption 
to the interfacial temperature of this quasi-steady con- 
duction-convection heat transfer problem, the matched 
CeMmpetbavuLrevatlongya: fin. on), =0..3 Te plottedian figure. (4:05) 
om aj log-log scale for Six values of, ¢ Gaal, Oye Wie Oe 
Pe Cpl Uma MOO Tt iShobvious, that winearity only 
occurs when the value of C is very small... In other words, 
a power law temperature distribution is a good assumption 
in the present problem only if the value of n is extremely 
small, that is, when it: is close to an isothermal fin problem. 
In general, the power law interfacial temperature does not 
exist for C equal to or greater than one. 
4,1.2-4 COMPARISON WITH EXPONENTIAL ASSUMPTION 

Lock [44] has shown that in two-dimensional free 
convection problems, the exponential surface temperature, 
=e for similarity, is merely an asymptote of the con- 
ventional power law distribution b=x™ when n approaches 
iifiniey ow Ln Order to checkiitheppossibmiityorL. the 
matched interfacial temperature following such an exponential 
form,: Logarithmic’ values of oe are plotted against x in 
PLGURewC4..0.).. 

Again, the linear relationship is achievable only 
when the value of C is much less than one. In addition, 
the profile for! C=l0sissalmost..a. straight, line near pihe 
HOOLTOLMEnNec tin. (x=1) >, this, 1ndacatesmienatswhen, the 


magnitude of C is much greater than one, an exponential 
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Figure 4.5 


C=0.01 


Quasi-Steady Interfacial Temperature 
Profiles for Checking the Possibility 
of Power Law Assumption 
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Figure 4.6 


C=0:50d- 


Quasi-Steady Interfacial Temperature 
Profiles for Checking the Possibility 
of Exponential Assumption 
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assumption may be a good approximation to the temperature 


distribution near the root. 
fe. S SOLUTIONS OF. OVERALL PROBLEM 
As far as the overall fin problem is concerned, 


PULereStcentres. on the total heat transter rate, the Nusselt 


number, and the fin effectiveness. 
4.1.3-1 TOTAL HEAT TRANSFER RATE 
if the base|area of fin is y2WD (See Ligure 2.2), the 


total quasi-steady heat transfer rate, Oar through thas (fpneas 


3 Paneelinatany 
Q. = BND ea (A2) 


and the heat transfer rate per unit area is 


ats (1) K.@, 465 (y) 
St oa Sar aa 2s Tite FRE RC EIS. ae 


therefore when both the dimensions and the temperature 

are fixed at the fin base, the total heat transfer rate 

of this quasi-steady problem can be expressed by the tem- 
d oS (y) 

perature gradient ace at. the tanwnoot. 

Reo LoOt Ob athas neat weranslLererabem Versus pDabamerer 
GMisvpresented insfigune (4.0) tore rave values of, y(0 7; 
O22) 0755 1 f0rdt ana One) ens -expected athe stronger che 
coupling of the convection and conduction problems, the 
higher the total heat transfer rate. A much more interesting 
result is shown in figure (4.8) where plotting total. heat 
transfer rate against parameter y on a log-log scale reveals 


aos (y) 
that linearity between eters a and ln y can be obtained for 
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Figure 4.8 
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every Crvalues* -Iniother words, the totalsheat transter rate 


can. be expressed by 


ag® (y) 
o = ay (4.4) 


WieresDOtie ra. and mM! are Lunctionseoreparamecer Cy / For 
SivectoalGie lanewin Ligure (4.0). im mists Slope and wat 
is the exponent of heat transfer rate at the intersection 
Of this line and the vertical line In y=0. Moreover, the 
Pana lel rcharacter of all) the strazghtlines in figure (4.8) 
suggests that "m" can be a common slope for all these lines 
and is hence independent of parameter C. 

In order to determine whether the slope mis a 
fixed constant for all values of parameter C within the 
range of this study, a least-squares method has been used 
to calculate the best fitted values of "a" and "m" for each 
straight line. Calculated results and associated fitting 


errors are listed in table below: 


Table 4.1 Curve Fitting for Total Heat Transfer Rate 


Versus y Plot with Different Slopes. 


C a(C) m(C) Sum of Least- 
Squares Errors 


0.01 0.00461 =Ie0a24 eer aley 4 
0.10 0.04518 ={eo324 One 
1.00 jesse? Lin OB Be don” 
3.00 0.90409 -1.0349 aoa 
5.00 1.29554 ese ey nes 
10.00 2.06882 Sees eis Was 
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these results suggest a much Simpler form of the 
empirical equation. The tabulated sums of least-squares 
errors are all in the order” of 1068 this indicates the 
goodness of these linear curve fittings on the log-log 
scale. In addition, all these slopes agree with each 
other to the second decimal place; that is, with less 
than 0.4% deviation. Based on this overwhelming obser- 
vation, a least-squares method was again used to 


calculate a common slope "m" and values of "a" for all 


the lines. The results are presented in Table (4.2): 


Table 4.2 Curve Fitting For Total Heat Transfer 
Rate Versus Parameter y Plot With 
Common Slope. 
mm = -1.03386 


€ a(C) | Sum of Least- 
squares, Errors 


0.01 0.00461 Teed ace 10 Ge 
0.10 0.04509 hse 
1.00 0.38095 DAs Out 
3.00 0.90538 1 Soe oi 
5.00 1.30000 Toss 
10.00 2.06170 Ometx to 


Sums of least-squares errors remaining in the order of 


Mino substantiates that an empirical equation takes the 
£orm aos (y) 


eye 


-1.03386 (4.5) 
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A further attempt was to express "a" in terms of 
powers of parameter C. A multiple linear regression 


technique [47] was applied to fit a polynomial in the 


form 


i Die) 7 es aie Pes ase A 5 
a Ch Dy = b,c + boC + 53C + b,c + bec 


The best fitted polynomial was found to be 


a(C) = 0.44066C - 0.067003C* + 0.007993C> = 0.00036377C" 


which yielded a total least-squares error of 0.015887 to 
the empirical equation (4.5). 

An alternative attempt was to fit an empirical 
Sxuacion in ap simprer) formsof 

da” (y) oF 

Ae SE aK: 

ae ace) Cuscyen) 

where function a'(C) was again expressed as a polynomial 
of parameter C. The result of the same multiple linear 


regression analysis program suggested a polynomial of 


a'(C) = 0.46123C - 0.070065c* + 0.00834c° - 0.00037781C* 


Based on this polynomial, the total least-squares error in 


fitting=equation (4.5a) was 0.589297 and the maximum 
percentage error was 10%. 
4.1.3-2 NUSSELT NUMBER 

For the classical laminar free convection problem 


adjacent to a vertical plate in a uniform gravitational 
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force field [39-43], conduction heat transfer in the 
solid wall or plate is not considered and the Nusselt 
number can be expressed by Nu = Nu(Gr, pr). 

Previous study [27,28] has Shown that when the free con- 
vection in the fluid is coupled with the conduction heat 
transfer along the wall and the gravitational force field 
pemre placed by ja centrifugal force (field suchwas\,in. the 
present quasi-steady problem, the Nusselt number is also 
aecuncrson Of parameter YY, and, Cra lf NG represents the 
local, Nusselt number at the fin root, thenfor Pr=0.72, 
Pewlses found that 


shes lie | 
Ni = Ra hii aE oi) (4.6) 


where Ra=GrPr, and Fi and F, are functions of the 
parameter C. 
For comparison, the local Nusselt number of 


present study is defined by 


: iia 
Sere tds) ls Silke [oo hab re : Ss 
Na = K, (ona Nu = E apy 7 (y,0) (4.7) 


Figure (4.9) shows that when 2O> (70) is plotted against 
parameter y on a linear scale, data points for every value 
of parameter C can be fitted accurately into least-squares 


line. In other words, 


-i/4 
ni SEE | = =o" (¥,0) = G, - Gy (4.8) 


where Gy and G. are functions of the parameter C. The 


calculated values of Gi and G, are lrsted in Table (4.3). 
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Figure 4.9 
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® Discrete Data 


Least-Squares Line 


C=10.0 


Nusselt Number of Quasi-Steady-State 
Solutions 
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Table 4.3 Curve Fitting For Nusselt Number Versus 


Parameter y Plot. 


cS Gy Go Sum of Least- 


0.01 fh ae Ge lay ex ome 
0.10 Oe5e6 ne Le 0.8 x 107° 
1.00 pecs 0.206 ioe TOE: 
3.00 0.1760 0.249 Dee oer Oue 
5.00 0.830 Ons oe atone 
10.00 0.929 0.297 OL Se Lome 


In order to make comparisons between equation (4.6), 


obtained from previous study, and equation (4.8) of the 
present investigation, Pr=0.72 can be substituted into 


the latter to yield 


ng = Ra/"(F! - FLY) Lye 


Here, Fs is the slope of the least-squares fitted line of 


the ay eae eeu VeDLot, and Fi Tswnociangrout che 


predicted value of nS/Ral/* Lor y=0sr 1n tLact, the static 
PeNVeOYrOo Lem .(35.)). 


Calculated values of Fi and Fs are tabulated in 


Table (4.4) together with their comparative quantities Fi 
and F, of previous work [27,28]: Apparently, there are 


some discrepancies between slopes EF. and Fy of the least- 


squares fitted lines. However, both Fi of the present 


od 


study and Pi of the previous work predict almost the same 


Ss 
value of Na/Rald * for the static fin problem. 


Table 4.4 Comparison of Least-Squares Fits 


Parameter Previous Work Present Study 


iS Fo F. Fi Fs 
On.01 Onsod OvnlekG 07,990 Oe 5 
Ono 0.401 07. 126 Ono L Ome 29 
P00 0.484 Oris 0.478 0.154 
3.00 O25. 73 Or 20'5 Ose IS On 266 
SAU) Orr oro 0.214 Or o2L OE2Z06 
10.00 0.696 OZ 28 Ort6.9°5 ORe272 


4.,1.3-3 Fin Effectiveness ¢€ 

Lock and Gunn [35] have defined fin effectiveness 
for a variable heat transfer coefficient h along the fin 
surface to be 

1 
fon (X) 0. (XK) ax 
= 42k 

A h(L) 6, (L)L Moot 

Since the numerator represents the total heat transfer 


Rate across, the thermal, boundarys layer iat one) side vot the 


fin, it must. be equal. to one half of the total heat transfer 


rate across the fin base, that is 
do (Lh) 
s dx 


L 4 x= cm 
fob (X)aX = -WK Cad) 


where W is the half-width of the fin base. Moreover, the 


denominator in equation (4:10) represents the local 
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convective heat transfer rate at the fin root multiplied 


Dy¥etnewcin length; hence 


h(L)0, (L)L = -LK, 0. (0,0) (Aga 2.) 


By combining equations (4.10)-(4.12) and also by 
using normalized variables (2.23) and transformations (3.5), 


the fin effectiveness becomes 


NS) 
_ yf 4(1+pr) pc 4 
Ci ¢| 4G Sco ete ee: 


In order to compare this with the fin effectiveness 
obtained in the previous study [27,28], Pr=0.72 is again 
used in above equation. For the sake of clearness, only 
dara LOtmsy=0. 1, 0.3 and O25 fare compared in figure (4.10). 
The present data, represented by discrete points, are 
in excellent agreement with the curves of previous 
results. 

4,2 PERTURBATION APPROXIMATION 

When the Froude number is not extremely large, the 
governing equations (3.17)-(3.19) and boundary conditions 
(3.20)-(3.21) indicate that solutions of this periodic 
problem may be affected by parameters Fr, Ki Ko, Cyorand 
y. By means of the series expansion (3.22), Fr is removed 
from all the perturbation equations and their boundary 
conditions. Nevertheless, numerical work in solving 
equations) )(5.23)-(3.28)) for the remavning four, parameters 


atavariols tLevels as Stil. expensive, only. selected 


P aa y 
pre F ; 
tag ay a 
, j, oN iz a ter 7 
> rs yy 
a t 
: i 
es 7 | o 
an ae 
~ ; Vie & 0 ¥ 
Cal r, 
i i 
i aedy 4 
+ ] - — 
mm! Apres a ies i 
in aie 
¥ y 
ie f e 
? a ' 
rae . 
F io 
7 i ais (os } 
5 roe 
a) ae See oe 
+ 
= 
4. 
i 
i 
‘ 4 : 
i a 
q 
= 
‘ 
F 
Sul 
i 
' St RR wore i) 
m Lin , U ao eens : 
i - oe yee Tae ca te , 
ao ale iy ey a ‘ ; 
4 ‘ MAES oY me ao, Ge 
. ~ Rad 5 i 
hae war fit ie i 
{ 
A iy 
it) = 
) [? a 
Ks f 
j ay + 
— =i) i 
Z | =: 
4 ‘ ie ‘is 
nh ; | 
v Lo } 
Pode 4 
; 5 Po . 
s ; j itl ' ] 
Cs ae , y : 
se iibs ee oe lag oe i fa te cl we 4 wey 2% a 7 id 
yy os ke, ly 2 ee in pn pers Wh a De —— J >" 
; salt i'n earls i : h * 
cj Y é ara, ‘ 
DI i] 1 a . 
7 
4 - 
i t 
: A, ' 
. 1 
eh 
i ri Ve s 
i ey, ‘ 
ues is 
pf , aa BA 
yy 
if ; 
- ' 


vo 
1 


: ge : } AS es m0 a a : re 7 Y P vi f° ‘ = in 
Day ; ; Te wr mn } gles 4 ad \ a a): , ’ 
Bal +.) GP A, A Pe! Ga e: 

F pene E/T 9 i ea, ay 


y i | 24 ine Ay ae ny 


€ 


FIN EFFECTIVENESS 


Zs 0 ® Discrete Data of the 


Present Study 


Least-Squares Line of 
the Previous Work [27] 


Figure 4.10 Fin Effectiveness of Quasi-Steady Solutions 
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parameter values are used here for economical reasons. 

EnNeerace OL SOluULTOn matching eatythe  solad—£luid 
interface was investigated for one Y value. From the 
quasi-steady results, the effect of y on the 3 Saeed order 
perturbation solutions was found to be considerably smaller 
than that of C. Basically, one representative level of 
y=0.1 was used throughout the investigation of this pertur- 
bation approximation, and only in some special cases, were 
solutions for y=0.5 generated. Similar to the preceeding 
study, the effect of C was reviewed at six parameter levels 
MOO eet ele e054 0 pang (LO0y meOr K,=1.0 and Kj=1.0; 
no difficulty was encountered during the integration of 
equations (3.23)-(3.28). Perturbed solutions of the 
conduction problem generally matched with that of the free 
convection problem along the solid-fluid interface within 
8-10 iterations. An exception was the case of C=10 which 
took 25 iterative loops to procure the convergence in two 
interfacial temperature profiles Bo and 20. On average, 
one iteration required approximately one minute computation 
on an IBM 360/67 computer. 

By setting C=1.0 and K,=1.0, the parameter effect 
OL Ky on the convergent rate of interfacial solutions was 
investigated at four levels (0.27 ,4.077 107 cand 100). 


Except when K,=100, solutions in general took 7-9 iterations, 


ie.) /—Jemanutes, to converge... “Anjsattempe for K,=100 was 


ee 


nee 
a we 


eiatal a a ie 
aa at a a 


unsuccessful due to the divergence of the solutions. An 


alternative attempt using K,=50 led to the converged solu- 


a 


LIONS anten 20st terative loops. WA further try using K=50 
and resetting C=0.01, reduced the required number of itera- 
tions to 5. This observation again emphasized the anti- 
cipaceds Hace vehat;the stronger the coupling) of conduction 
and Eree convection problems, the more difficult the 
convergence of their solutions. 


The effect of parameter K, on the solution matching 


2 


at the interface was also studied. When C=1.0 and K,=1.0, 


1 
SsoLlucioms: converged. within 7 iteratronss forall four para-= 
meter levels of K, (Om Or Omran Cer EOOD ES ae LI taCtr ie 
further attempt for K,=1000 also produced converged result 
within 1l iterative loops. It was therefore concluded that 
when C=1.0 and K,=1.0, the convergent rate of the pertur- 
bation solutions was rather insensitive to the parameter K.. 
4.2.1 /SOLUTIONS (OF ‘BOUNDARY -LAYER PROBLEM 

Prior to investigating velocity and temperature 
profiles of the recombined (i.e. perturbed) results, it 


Ms uiogieal to explore some characteristics of, the pertur— 


Daetion solutions. 


The amplitudes of the first-order perturbation solu- 


ELOnS: ae 2, ie and $, for a typical example of y=0.1, 


C=0.01, K,=1.0, and K,=100 are presented in figure (4.11) 


80 


letes ey ey ie i 0) ab an i i ahs 
if we - a5 i} a Fe ¥) ad H y 
4 Meher metals rie ier fs > aera be fa Jad cee i 
+ 5 v Z c , r ir 
i 5 re by es ; 
ry 
Pi s © F 
J it 
Sy 
o- PM 1 
> i i . 
2 Ky 
}iee. 
. : ' iS 
ae ‘ be ie (ta 
2 - 1 
‘i 
‘ 
¢ : 
7 os 
ot Zé 
Si i 
i J te fe 
i } ty 
= 
i “<2 
b 
- & 
j t Mi 
is F a ' 
= i 
: & 
; 
af 
f 
/ 
st \ 
\ 
’ a 
y f want 
i 
ee 
J ‘ 
\ 
ps : 
i sh =F 
i Boas ; : 
‘ 
os ton 
= 7 
i 
ee } vs ~ 
s 3 if 7 
‘ ro A 
, ‘ , or ital Fie) Pow 
i \ i : J ix” wa ae P) Pe daey! wh id 
? ' i ; , Aw | 
; UL TMA ate Ae 
3 “4 ak 
i 
ft 
1 
f , ‘ " 
i 
=) a 
Ay si 
x ; Ps 
tg = = Gs Le fess) Le 
1 y, 
b - ; 
: i 
: i = r ao 
Ee Pa ey 7 Git eS a 8 ie 
& i i i } a i 
Tate? ; 
Re Le 
ETS 
” = 
E 
; in 
- cs 
} 4 


0.00 pm ee eee eee 40 
e Beate 
g> a 
7 
vA 
/ 
fe Velocity 
-0.02 / — ae Lemperature Rts 
-0% 04 SUPEES Crip Ecma whom x=04:2 0.24 


DaptoOir x=i70 


=O 6 O36 

f 

n 
—-0.08 0.08 
eA D 828) 0.00 
Syste. -0.08 
-0.14 —-0.16 

0.00 2.00 4.90 6.00 8.00 10.00 


Figure 4.11 Even Part of the Perturbation Solutions 
for the Boundary Layer Problem 


avd (aU2 i ethne odd) parts of) tnesperturbation to the steady 
solutions are approximately ten times smaller than the even 
parts. Examination of all the cases under this study shows 
that this observation is generally true except in the cases 
of K,>0f1], where both odd and even parts of the perturbation 
Solutions are of the same order of magnitude. In spite of 
these relative magnitudes, it is more important to determine 
whether they can introduce significant periodic effects 

with a reasonable degree of accuracy. 

In general, the truncation error of a first-order 
perturbation series expansion can be estimated by the 
Magnitude of the second-order terms, Accordingly the 
accuracy of the present results relies on the product of 
Pr? and the second-order perturbation solutions. On 
account of the fact that both steady and periodic parts 
of the solutions have been properly normalized in 
Appendix B (during the derivation of governing equations), 
the magnitude of the second-order perturbation results is 
Offorder one, ana vthe accuracy of ehis first-order pertur-— 
bation approximation: is therefore in the order of Pouce 
Slimane rancat.on, error of 0 [1077] is the acceptable limit, 
ten is then the smallest Fr value to be used in studying 


these perturbation solutions. 
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Prgqure: 4.12 ""Cdd "Part oP7themrerturpation‘Solugions 
for the Boundary Layer Problem 
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Figures (4.13)-(4.18) show some examples depicting 
the periodic features of boundary layer solutions for 
different parameter values. None of them shows a noticeable 
periodicity in the fluid température when Fr is equal to or 
Greater than 10; on the other hand, velocity profiles for 
Be 0vallv reveal some periodic feature: in the vicinity of 
their maximum values. Figure (4.18) shows that this 
periodicity disappears in both temperature and velocity 
profiles when Fr=100; it would appear that the problem has 
reached its quasi-steady state when Fr approaches this 
order of magnitude. 
4.2.1-1) EFFECT. OF) PARAMETER ¥ 

Boundary layer solutions at the roots of two 
different fins (y=0.1 and 0.5) are compared in figure 


a2) nom C=0.01, Ke=1..0, K,=1.0, andsrr= 0p weOtimathernan 


1 
and momentum boundary layers are seen to increase their 
thicknesses with parameter y; but as far as the normalized 
quantities are of concern, the greater the parameter y the 
smaller the magnitude of the maximum velocity. Although 
the non-periodic temperature across the boundary layer for 
V=055 usmoenerally higher thanythat: for y-0.8, the steeper 
temperature gradient for the latter indicates a higher 
local heat transfer rate at the fin root when the parameter 
y is smaller. 

Figure (4.13) also shows that solutions vary with 
themtinepos@fion. athe inclination or  herrotating (fin is 


indicated by the angle starting from 0° when the fin is 
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Bagure 4715 Perturbed Boundary, bayer Solutions for 


Fr=107.,C=2; K,=1, K,=1 and. y=0.1 
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Figure 4.16 Perturbed Boundary Layer Solutions for 
Fr=10, C=l1, K,=90, K,=1 and y=0.1 
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Figure 4.17 Perturbed Boundary Layer Solutions for 


Fr=10, C=l1, K,=1, K,=100 and y=0.1 
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pointing downward and rotating counter-clockwise. When the 
fin 1S making a 0°-angle with the vertical line, the buoyancy 
Serlects Or Gravity sand centritugaletorce concur; the resultant 
velocity fh is hence greater than that for the 180° DOSt ELON, 
in which these two driving forces oppose each other and their 
total buoyancy effect decreases accordingly. No difference 
PieticevelOci ty protiies for the finjjat the two horizontal 
positions has been noticed. 
4.2.1-2 EFFECT OF PARAMETER C 

Figure (4.14) discloses the effect of parameter C 
on the thermal and the momentum’boundary layers. The 
thicknesses of these boundary layers are insensitive to 
parameter C, although both the temperature and the velocity 
are greater when C is smaller. Moreover, the temperature 
gradients at the solid-fluid interface (n=0) suggest a 
Higher heat transfer rate from solid to £luid) at the £in 
root when the parameter C increases. This is logical 
because the more conductive the fluid (he. the greater C) 
the higher the heat transfer rate. Again, fluid velocity 
FOr Mcre vt 11 0° is greater than that at Oo} and no 
deviation is noticeable between the velocity profiles for 
the fin at 90° and 270° positions. 
4.2.1-3 EFFECT OF PARAMETER Ki 
By using the same parameters C, Koy andsEry athe 


velocity and temperature profiles of the fluid at x=1.0 


eae a ae 


ei . ait i 


ie 


bs 


Lor K,=1.0 and 50.0 are presented in figures (4.15) and 
(4.16) respectively. For K,=1.0, no difference is noted 
Mieciewve locrry profilie=when the £inwvis vat) two horizontal 
positions, (couyesponding (to. 90f,,and 2701 in) the drawing). 
Soiswagreement. in velocity profile is’ mainly” due to the 
extremely small value of ee in the perturbation solutions 
when K,=1.0. However when K,=50, figure (4.16) shows that 
the velocity fh for the fin at 90° deviates appreciably in 
its peak value from that at 270°. On the other Wana) secoues 
welocity profile hardly varies as the fin rotates from 0° 
to 90° or from 180° to 270°. This temporary steadiness in 
fluid velocity is because the odd part os of the perturbation 
solutions has the same order of magnitude as the even part 


~ 


aay and generates noticeable effect on the solution when 


the parameter Ky increases. 
The rate of solution convergence can be affected by 
parameter K,. Mathematically, this Ky 


intensity of coupling between the perturbed boundary layer 


value represents the 


equations while physically it indicates the significance 

of transient effects on the boundary layer solutions. As 
Mentioned berore, the difficulty in achieving the inter- 
facial convergence of the perturbed solutions increases 

with parameter Ki, but decreases more rapidly with parameter 
C. In other words, the parameter C plays a much more 
important role than Ki in the convergence of interfacial 


solutions; however, a large K, -value does aggravate the 


convergence problem when the parameter C is not much smaller 


Chem sumdccy. 


4.2.1-4 EFFECT OF PARAMETER K, 


Tikesrtigure.(4915)o for Kj=1, figure (4.17) depicts 


Thee per curbedsboundary* layer? solutions tor KA=100 while the 


2 
other parameters remain unchanged. Almost identical 
solutions from the above two figures reveal that this 
parameter K, has no influence at all on the perturbed 
boundary layer solutions. 
4.,2.1-5 EFFECT OF FROUDE NUMBER (Fr) 

The larger the Froude number the smaller its effect 


on the solutions. Figure (4.18) gives the boundary layer 


solutions at x=1.0 for Fr=100. Comparison of these solutions 


with those presented in figure (4.15) discloses that the 
periodic feature vanishes, as expected, when the Froude 
number becomes very large. In fact it is reasonable to 
predictothat -thercsolutronian figure (4.18) is very close 
to the quasi-steady results for extremely large Froude 
number. This prediction is proved to be correct as 
solutions in figure (4.18) match completely with the 
quasi-steady solutions for y=0.1 presented in figure (4.2). 
442-2 o sSOLUT TONS OP. CONDUCTION PROBLEM 

Parameter (C’has va «significant effect on the conver- 
gence of the perturbation temperatures along solid fin. 
The magnitudes of all the perturbation temperatures o, and 


oo obtained in this study are so small that after being 


gS 


Ne 


ids 
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divided by the minimal Froude number of 10, none of them 
can give rise to any noticeable change in the quasi-steady 
solution he along the fin. Therefore, none of these 
parameters noticeably affect the solution ube in the conduc- 
tive fin problem. Nevertheless, experience shows that it 
PomMOuesiite1Oultto Saristyy thevconvergentecuLterion of 


107° 


for these perturbation temperatures o and $C when 
the parameter C becomes larger. This implies that at 
least one of the above perturbation temperature either 
increases its magnitude or change profile with C. 
Convergent rate of the perturbation temperature 
along the fin decreases inversely with magnitude. Figure 


~ 


(4.19) depicts solutions as and “Bs POrsy=0 . pe Keel at 


if 
three levels of parameter C: 0.01, 1.0, and 10.: When C=0.01, 
both o. and o are very small and their magnitudes decrease 
monotonically from fin tip (x=0) to fin root (x=1); these 

two favourable features contribute to the rapid convergence 
Stmchenscolutcions. at, the isolid-filuid interface. "As the 
parameter C changes from 0.01 to 1.0, the magnitudes of o 
and o. increase approximately fifty times without losing 

the aforesaid monotonically decreasing feature; this 

increase in magnitude results in the requirement of more 
iterations for the interfacial convergence. 


Shapes of the perturbation temperatures ae and pe 


Significantly affect the rate of solution convergence. 
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First-Order Perturbation Temperatures 
Along Soliad’rin 


an obvious change in the profiles of these perturbation 
temperatures appears when C=10. In spite of keeping the 
Same order of magnitude as when C=1, both = and $ 
increase with distance x from the fin tip to their maximal 
magnitudes at certain (but different) x-values, and then 


decrease inversely with x to zero at the fin root. 


phuS loss) of monoconically warying (characteristic in ithe 
magnitudes of 2 and 3 is possibly the major hindrance 
to the rapid convergence of these perturbation solutions 
at the solid-fluid interface when the parameter C is much 
greater than unity. 
4.2.3 SOLUTIONS OF OVERALL PROBLEM 

It has been found that when the Froude number is 
not less than ten, the temperature of neither fluid ee 
solid fin experiences any noticeable disturbance from the 
perturbation solutions. An alternative investigation may 
be undertaken on the possible disturbance in temperature 
gradient and heat flux which, as poner geeate in section 
Pes are, close ty related tothe totalwecatstransfer rate, 
the Nusselt number, and the fin effectiveness. 
4.2.3-1 TOTAL HEAT TRANSFER RATE 

By referring tovequation (4.3); the jheat transfer 


rate per unit base area of the fin is 


K_9O 9, Cy, 7) 
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That is, the total heat transfer rate across the fin base 
at time < is proportional to the temperature gradient 


JO Cy, -1) 
fe) : 
amet! which can be further approximated by the pertur- 


Daclon= result: 


eo u ao” (y) aes dd. (y) ad (y) 
O& | SIS Aa aes [—ge— cos (2nT) + —g¢— sin(2mt)] 


(4.14) 


Figure (4.20) demonstrates the periodic variation 
of the root heat transfer rate as the fin changes its 
position angle a while rotating (the fin is projecting 
downward when t=0 or 1). It is clear that if Fr increases 
to 100, the total heat transfer rate approaches that of 
the quasi-steady solution. When Fr decreases to 10, the 
Eluctuating amplitude of this total heat transfer rate 
increases to approximately 2% of the quasi-steady solution. 
A dashed line for Fr=2 is also presented in figure (4.20) 
though the accuracy of this first-order perturbation 
approximation is considered to be extremely questionable 
whens tis Ft value’ is Jess than’ 100. the plot, however, 
still indicates that a fluctuation in the order of 10% 
is plausible for Fr=2% 

All the curves in figure (4.20) have some common 
features which may be used to determine the proper initial 
time level for the direct integration approach. First, 


every curve is fairly linear when t varies from 1/8 to 3/8 
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and PEOmns 7791 tO. 7/8. In addition theroeak point and 
valley wpoumteor every curverare always situated at t=0 
and 1/2, respectively, and all the curves intersect the 
(Udsi=oceady result in the) vicinity onmtwo pornts< t=1/4 
andnsfjge Aliens “suggests ‘that either t=1/4 or 374 may 
be chosen as the initial time level for the direct 


integration of governing equations (3.17)-(3.19). 


4.2.3-2 . NUSSELT NUMBER 
From equation (4.7), the Nusselt number at the root 
of a rotating fin for this periodic heat transfer problem 


can be deduced: that is 


4(1+Pr) 


where o. can be approximated by 


i/ 4 
Nu = =| ate 0, (y00,7) (4.15) 


o (y0rt) = O(y,0) + Feld, (1,0) cos (2mt) + 8, (y,0)sin (271) )} 
(4.16) 


A plot of this temperature gradient versus time Tt 
(and position angle a) for several magnitudes of Froude 
number is presented in figure (4.21). By comparison with 
Prgure (4, 20), thiss plotireveals thatthe Nusselt number 
at the fin root has a lower percentage fluctuation than 
the total heat transfer rate in the present case. Also, 
aki the points: Of sinterest, oye. schnempeacmpolncs, the 
valley points, and the points where curves intersect the 
horizontal line for the quasi-steady solution, seem to 


ehectectiaqnely, rightward |(102)) 1necigureg(4..21)\7 9 (inwspite 
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of thas, the nearly linear feature’ of all the curves when 
the SrOLatingmran Changes Lts positions £com) tal/4 ito 7=3/8 
and from 1=3/4 to t=7/8, again suggests that both 7=1/4 
and t=3/4 may be the proper initial time levels for the 


eLEccupincediatl1On, approach. 


iMgcact, /evenuitathis ploeror t(y,.0), 7) versus T 

oo ott) i 

does NOG Support that of ‘Kone eee in determining 

the initial time level, the recommendation from the latter 

305 (YT) 

should be the proper one because —4F represents a 

physical measure of the overall problem while fy, 0) only 

provides a local indication. 

202.5—5 FIN EFFECTIVENESS 


A significant amount of phase-shift appears in 


figure (4.22) where the fin effectiveness ¢, defined by 


Ses 


G 3pr aay me oe Ma 28) (ah lg 


is plotted versus time T and position angle a of the 


¥ Nn i 


rotating fin. Apparently,-all the points of interest show 
the same excursion to the left of the graph. For example, 


the fin effectiveness predicted by this perturbation 


solution reaches its maximum value approximately 30° before 


the fin rotates to its vertically downward position (t=0 


or 1) where the maximum total heat transfer rate takes place. 


30 (Y,T) 
Unlike a and De alas which are real 


entities directly involved with the integration of the 


Governing, equations, fin eftectiveness: defined by equation 
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(4.17) is merely an intuitive appraisal of the problem; 
neither physically nor mathematically can it affect the 
cicece PNLegration-towany oxten tan Bom chisereason use. of 
figure (4.22) is not recommended in the determination of 
initial time level for the direct integration approach 
though both t=1/4 and T=3/4 are still indicated by this 
plot to be two time levels of satisfactorily small curvature. 

Again, the fin effectiveness for Fr=100 fluctuates 
almost imperceptibly about the quasi-steady result in 
figure (4.22). This fluctuation increases to about 1% and 
3.5% as the Fr value decreases to 10 and 2, respectively. 
Rie oChCdy, «SOLUTIONS 

Pn mexact = “Ssolbutiron) ine thaissctudy) refers: to the 
numerical solution obtained from the direct integration 
of equations; (3.17) -(3.19).,..General.guidelines.for this 
Girectswintegration, have) been described) inisectioni3).)2)2 
through some examples depicted by figures (3.1)-(3.6). 

The suggested use of first-order perturbation to 
approximate the exact solutions at a proper initial time 
level in this work postulates that the variation of 
Froude number would not cause significant phase-shift in 
the solutions; that is, cases similar to either figure (3.1) 
or, (3)..3)\ are assumed...» This assumption) is) later verified 


a posteriori. 
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Because of the expensive numerical work expected 
PnetiiSsrappLroacn Of ‘direct inbecraciens the study centred 
on the effect of Froude number, the most important para- 
meter of this periodic problem. Exact solutions were 
obtained at four levels of Froude number (2, 10, 100 and 
10>) for a representative set of parameters: veal Pra 72, 
C=0.01, K,=1.0, and K,=100. 

In order to reduce the computing time, integrations 
were carried out by using’ the step sizes of An=0-1, Ag=0.01, 
and At=1/8. By virtue of the first two step sizes, quasi- 
steady and perturbation solutions could be secured within 
seven and three minutes of computation on an IBM 3640/67 
machine. 

In the direct integration of the governing equations 
(3.17)-(3.19) for a large Froude number, the proper initial 
time level was determined from the perturbation result. 

In this study, the predominance of the even part of the 
perturbation solutions over the odd part hinted at the 
possible resemblance between the exact solutions anda 
cosine curve. Moreover, the assumption of no significant 
phase-shift among solutions for different Froude numbers 
suggested approximately linear solutions with respect to 
timein the vicinity of t=1/4, which was “thererore "chosen 


as the initial time level for large Froude numbers. 
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For comparison, quasi-steady and perturbation results 
Wevemseparate.y, used asm the initraresomutronvatet=l74.~ By 
means of the quasi-steady result, timewise convergence was 
achieved for acl arter 135 minutes erscomputation, but 
failed completely for Fr=100. When the perturbation 
approximation was used, however, solutions in both cases 
(Fr=100 and 10?) were obtained with 25 minutes of computing 
time. Therefore, a search for the perturbation result 
prior to the exact solution is both necessary and economical 
even 2f the Fr ovalue “is not small. 

The initial time level for all the small Froude 
number cases (Fr=2 and 10) was also set at t=1/2. In the 
present study, the exact solutions for large Froude numbers 
all deviated from their perturbation approximations in the 
manner demonstrated by figure (3.1); namely, no appreciable 
phase difference was revealed. It was assumed that same 
phase behaviors held between the exact and perturbed 
.solutions for a small Froude number. | 

Numerical experiments were conducted to choose a 
proper initial solutions for the small Froude number. 

With the exact. solutions) tor, Fr= 100 as the initial solutions, 
the periodic solutions were matched to better than the third 
decimal place within 40 minutes of timewise iteration. The 
overall computing time for a small Froude number was in 


the order of 75 minutes. Both quasi-steady and perturbation 
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results were also used separately as the trial initial 
SsOoluticonsein this direct integration approach tor small 
Froude numbers. The former often led to solution diver- 
gence or numerical overflow while the latter produced 
convergence after 200 minutes of extremely long calculation. 
Undoubtedly, direct integration of the governing equations 
for a large Froude number before solving the case for 

small Froude number, would shorten a prodigious amount of 
computing time. 

The demonstrated approach for the direct integra- 
tion of equations (3.17)-(3.19) assures fast timewise 
convergence of the periodic solutions. Notwithstanding 
all endeavors and successes in saving computing time, the 
minimum requirement of 75 minutes for solving a small 
Froude number case may still make the whole scheme seem 
economically unattractive. Closer examination reveals 
that this long computation Oey reflects the complexity 
of the coupled heat transfer model at each time step. 

In fact, periodicity of the solutions can be obtained 
after three timewise iterations. 
4;3.1 SOLUTIONS OF BOUNDARY LAYER PROBLEM 

To compare the exact solutions of the boundary 
layer problem with their perturbation approximations, two 
sample solutions for Fr=10 are presented in figure (4.23) 


and, 424). “Apparently, ‘there is iveny little difference 
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Figure 4.23 Periodic Boundary Layer Solutions Obtained 
from Perturbation Approximation for Fr=10 
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Figure 4.24 
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Periodic Boundary Layer Solutions Obtained 
by “ther Diresctmintegrartion Method for Rr-10 
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in the temperature profiles 6 between these two graphs. 
In contrast to this, the velocity approximated by the 
perturbation results in figure (4.23) is generally higher 
eiengetoeexact Value in figures (4-24 )edespitesthe fact 
that they both fluctuate in nearly the same magnitude 
and reveal the same boundary layer thickness. As far as 
the accuracy of the perturbation approximation is concerned, 
however, the greatest velocity difference of three per cent 
between these two sample solutions is still within the 
truncation error, Oomom. of the perturbation results. 
Figure (4.24) depicts physically sound velocity 
protiles my Aswvthe, fin is projecting downward, ive. Om 
concurrence of the gravitational and the centrifugal forces 
induces a/£luid velocity higher than that of iene where 
Opposition of forces reduces the net buoyancy effect on 


the fluid. No appreciable difference exists between 


the velocity solutions at the two horizontal positions 


(90° yang? 2,707) 
The fluctuation in fluid temperature profile is 
very small. For this reason, only the limiting temperature 


profiles: at two vertical win positions: of 0° and 180° are 

plotted#ain figure (4.24); temperature profiles. ati other 

fin positions are virtually enclosed by these two limits. 
The boundary layer solutions obtained by the direct 


integration method for Fr=2 is shown in figure (4.25). 
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Figure 4.25 Periodic Boundary Layer Solutions Obtained 
by the Direct Integration Method for Fr=2 
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Thestemperature profiles vary slightivy with the position 

of the rotating fin but the velocity profiles fluctuate 

so intensively that they even present noticeable differences 
Sotne two ener, Zonta l fin (posit ONS ser ineparti cular, -the 
Momentum boundary layer at the position angle of Det) ae 
thicker than that of 90°. The solutions also indicate 

that although the maximum velocity varies significantly 
when the fin rotates from 0° to 90° (or from 180° to 270°), 
the thickness of this momentum boundary layer stays nearly 


unchanged. 
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4.3.2 SOLUTIONS OF CONDUCTION PROBLEM 

Figure (4.26) shows the periodic fluctuation in the 
solid fin temperature obtained/by direct integration of 
eclations @at7)— (3.19) forwr=2.0 9 When the Heated rotating 
fin is projecting downward ae concurrence of centrifugal 
and gravitational forces enhances the cooling effect of the 
ambient fluid and, as shown, results in a lower fin tempera- 
ture profile. On’ the other hand when’ the rotating fin is 
projecting upward amore opposition of the two body forces 
reduces the cooling effect and accordingly leads to a higher 
fin temperature distribution. Similar graphs have been 
constructed for other Fr values and, as expected, the 
difference between these temperatures at two vertical fin 
positions varies inversely with Fr and eventually disappears 
when Fr is greater than 100. It is also shown in figure 
(4.26) that when the fin is in the two horizontal positions 


(90 


and 2708), the difference in the corresponding fin 
temperature is negligibly small. 
433..3)) SOLUTIONS, OF OVERALI: PROBLEM 

Periodic variations in the total heat transfer 
rate, the Nusselt number, and the fin effectiveness for 
exact solutions are presented in figure (4.27), (4.28), 


and (4.29) attour levels of the Broude number, (2, 10,7 100 


and 10°) when C=0.01, y=0.1, Pr=0.72, K,=1.0, and K,=100. 
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Figure 4.26 


Periodic Fin Temperature Obtained by the 
Direct Integration Method for Fr=2 
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For comparison, the perturbation results are also plotted 
as sOdsnecm tunes fOr oOtNerrT=2 mandeslOr 
4.3.3-1 TOTAL HEAT TRANSFER RATE 

Brommequations (442 )mandu (4s) eetnestotal heat 


Evanster srate © for this peniodiclproblem can be derived 


as K_0 Ou (¥ 5) | 
Q(t) = | 20 > | co Aides See Care 82) 


R, aE 


As all the bracketed quantities in the above equation are 


Sect COns tant, for a given problem, this total imeatsetranster 
ae. (ae 3%) 


o€ 


rate can be expressed by the temperature gradient 
ee sene = TOocLof the solid. fin. 

Several characteristics of the total mys transfer 
rate for a sample problem are disclosed in figure (4.27) 
where the temperature gradient is plotted against time T 
and position angle a of the rotating fin. In agreement with 
the perturbation approximations, maximum heat transfer rates 
(eak poin ts)) of -exact. solutions for ir=100).70.and #2, all 


(t=0 or 1) while their minimum heat transfer 


occur at a=0 
rates always take place at a=180° (t=1/2). Moreover, all 
the curves intersect the quasi-steady result in the vicinity 
of a=90° (t=1/4) and a=270° (t=3/4). The nearly linear 
sections between tT=1/4 and 3/8 of all the curves shown in 
figure (4.27) verify that the use of t=1/4 as an initial 


time level in obtaining exact solutions for these cases is 


perfectly acceptable. 
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It is also worth noting that magnitude of the 
Froude number can affect the total mean heat transfer rate 
of the problem. When the Froude number is much greater 
than one, the total heat transfer rates indicated by the 
exact solutions in figure (4.27) fluctuate almost symme- 
trically about the quasi-steady result; namely, their total 
mean heat transfer rate equals the quasi-steady solution 
to a close degree. However when the Froude number is of 
order one, this mean value becomes significantly less than 
the quasi-steady result. 

Discripancy in the total heat transfer rate between 
the exact and the perturbed solutions increases inversely 
with the Froude number. In general, a decrease of the 
Froude number undulates both the perturbation result and the 
exact -soOLuLTOn ; DUL to a greater extent in the latter case. 
Figure (4.27) reveals that the maximum discrepancies between 
the exact solutions and their first-order perturbation 
approximations are respectively 0.5% and 7% for Fr=10 and 2. 
In comparison with the estimated corresponding truncation 
Sig eere yy Eee o[fr “I, Of this first-order, perturbation: series 
expansion, these low discrepancies suggest that the total 
magnitude of all the high-order perturbation terms is rela- 
tively. smallain tht Ss study <andmtiemuser Of tne mii ns t-Order 
perturbation approximation alone is acceptably accurate in 


most cases provided the Froude number is not too small. 


4.3.3-2 NUSSELT NUMBER 

Equation. (4.15) indicates ¢that ton a given,problem, 
the Rayleigh number (Ra) and the Prandtl number (Pr) are 
fuxed and ythepNussel.t number atsitheuwoot Of .a rotating, fin 
can be expressed by the temperature gradient o(y,0,T) of 
Eheyriuidearmthe soolid—-tluidgintertace.  wihis temperature 
gradient is plotted against time t and position angle oa 
mnetegure (4.28). 

Several characteristics of the time-dependent 
Nusselt number are observed. In comparison with figure 
(4027), atl, thexcurves iin sfigure:.(4 428) stake an, excursion 
to the right. Nevertheless, linearity still appears between 
T=1/4 and 3/8 for every curve; this trend again substantiates 
the wisdom of choosing t=1/4 as the initial time level in 
the direct integration approach. Similar to the total heat 
transfer rate, the Nusselt number undulates its fluctuation 
but decreases in its mean ite as the Froude number becomes 
smaller. 

It is found that the verturbation approximation is 
less reliable when Fr=0O[1]. As far as the magnitude of 
fluctuation is concerned, the perturbation results of both 
figure (4.27) and (4.28) estimate a lower percentage change 
in. the Nusselt number than in the total heat transfer rate 
for Fr=2; exact solutions, however, manifest an equal 


percentage: variationyin bothventitiesmand hence, disprove’ the 


observation. deduced: from perturbation results. iin figure 
(4.28), the perturbation results for Fr=10 and 2 can deviate 
up to 1.5% and 11% respectively from their approximated 
exact solutions. These deviations still concur with the 


estimated truncation error, O[Fr7- 


J, of the perturbation 
approximation. Apparently, the use of the first-order 
perturbation series expansion to approximate the exact 
Solutronsiiat therchosen initial’ timeslevel’ t=1/4 in this 
study is sound if the Froude number is much greater than one. 
4293..3-3 PEIN EFFECTIVENESS 

Povacions (4017) kofesection 4). 362-3016 Usedsin 
determining the fin effectiveness ¢€ of every directly 
integrated exact solution. The calculated results and their 
approximated values for both Fr=2 and 10 are plotted versus 
dimensionless time Tt and fin position angle a in figure (4.29) 

Several distinct features are present in figure (4.29) 
First, none of the periodic fin effectivenesses resembles a 
cosine curve; all the peak and valley points deviate from 0° 
(t=0 and 1) and 180° (t=1/2) to an appreciable degree. For 
Fr=10, the exact solution indicates a fin effectiveness 
fluctuating about the quasi-steady result (Fr=107>) rea ee 
to a negative sine curve, i.e. with the valley at 90° and 
peak at 270°. In other words, the phase angles between a peak 
and ahfollowing valley; or vice wersay is always *haliga perrod 


(180°); this also holds for the perturbation results for both 
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Pr-2 and 10. The exact solutions for Pr=2, however, present 
completely different features for fin effectiveness, which 
takes about one third of a period to rise from its: valley 
to its peak but needs the remaining two thirds of the period 
to subside to the valley point of the next cycle. In addition, 
this fin effectiveness of the exact solutions for Fr=2 reveals 
a mean value higher than the quasi-steady result although its 
perturbation approximation predicts a lower mean value than 
the latter. 

Despite all the differences in phase angle, shape, 
fluctuating magnitude, as well as in the mean value, maxi- 
mum descrepancies in fin effectiveness between perturbation 
result and their corresponding exact solution are 0.8% for 
Pi LO and 66 Lor Fr=z2. 
4.4 _GIMITS OF, APPLICABILITY 

ROD Teabisityveor. the present result is subject Oo 

some basic constraints of this problem. Assumptions utilized 
in appendices A and B to derive the laminar boundary layer 
equations (B.34)-(B.37) impose the following limits on the 
solutions: 

(1) Limit for laminar free convection phenomena. 

In order to avoid the development of turbulent flow in 
conventional free convection problems, i.e. gravity is the 
only body force for the entire system, the Rayleigh number 
should be less than 10%. By assuming that this limit on 


the Rayleigh number also holds for a laminar free convection 
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(2) 


(4) 


problem in a centrifugal force field, the Rayleigh number 


defined by BRw ‘0 L* 
aN, : 
should satisfy the following: 
Ray< 0. a (4.19) 
Limit for boundary layer flow. 
A boundary layer flow will take place when 
Rae one (4520) 


Limiter or 1 nStgni fli Cant: Corio lisetorce: 


The Coriolis force effect in the momentum equation 


(B.11) can be ignored if 


Ek7Ra73/4 aan Copia 


where the Ekman number is defined as 
yay) 
oe ear) 
For convenience, the above condition is replaced by 


Eko Rags hoe ee (AR O18) 
Limit for negligible viscous dissipation. 

All the viscous dissipation effects, represented 
by the last three terms on the right-hand-side of the 
energy equation (B.13), become vanishingly small if the 
ostrach number defined by 
pRo*L 
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(oy)? Dam t tor novirevers al flow. 
When the centrifugal acceleration at the fin root is 
gneater ethan’ the qravetationalwacceleration, i.e. 
(R= )G) eee (45523) 
the fluid in the boundary layer will always flows 
radially inward; no reversal flow will take place. 
(6) Limit for Oberbeck-Boussinesgq approximation. 
Applicability of the Oberbeck-Boussinesq approximation: 
OP =FipR el Bo)a ae, (ASS?) 
relies on the assumption 
BO < can lime, (A.6) 
Or more specifically, 
BOM Orel (4.24) 
Actual constraints on the applicability of the 
solutions given here can be clearly illustrated by delineating 
these limits on a properly chosen coordinate system. If 
air at a temperature of 68°F is the ambient fluid of a 


6-inch rotating fin,and 0.=10°F, the estimated value of 


10 yf 
BOT aeAGORGGe iA 
automatically satisfies the inequality (4.24). Furthermore, 


substitution of the physical quantities [48] of air: 
Ke = 236 x NOGiet tay secur 
Vi 160 ex LlOm sete sacs 


and Co = 6000 ft2/sec?-°R 
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into the remaining inequalities (4.19)-(4.23) yields the 


following more rigorous constraints: 


Ro < 12595 (4.19a) 

Rusa 20. Lo ee (4.20a) 

Rw2/3 > 6.232 - (Ate? ey) 

Rw? < 633600 (ARG20 a) 
and 

(R-0.5)w? > 32.2 (Ae2 34) 


These inequality Re nae are shown in figure (4.30) 
where 1lnR is plotted versus Ilnw. The shaded side of each 
constraint contour indicates the region where the solutions 
illustrated can not be applied. Applicability of: these 
solutions /is. limited to the, dotted region in figure (4.30). 
The actual constraints are, therefore, determined by limits 
Fs) ana (5); ise Muinequall Eves (4.119 )) Ne 4.2%) 0 and (426)) 
in figure (4.30):7 ditiereng portions \of the dotted 
region corresponds to different possible applications of the 
presented results. Solutions for the upper portion of the 
dotted region can be used in a re-entry environment because 
of large radius R and low angular velocity w. On the other 
hand, solutions for the dower" postion, which ws ichagacterizeca 
by small radius and medium angular velocity, may only be 
applied to a centrifugal machine rotating at relatively low 


Eome s bt isisnown in figure (4230) Bthat the highest) angular 
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Figure 4. 30 
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Rw*>0.16 
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velocity, for the dottedyregionyshould be less) than the angular 
velocity associated with the intersection of limits (1) and 
(5). By solving the simultaneous equations: 

Ree 1556 
and (RA04.5) cepi=93 27,2 Sy, 
which are deduced from these two limits, the eee velocity 
for the intersection is found to be. 

ol = 56 rad/sec = 535 rpm; 
such a low limit on the angular velocity indicates that 
solutions illustrated have no application to high speed 
turbines. In addition, based on the inequality (4.19a), 


Enhes solutions can not be applied) tovasrotatingusystem 


Wath the centrifugal acceleration greater than 50g. 


CHART E Ravi 


CONCLUSIONS AND RECOMMENDATIONS 


Sele CONCLUSIONS 

The two chief objectives of this work have been 
esuccessiully achieved. First, a periodic, heat eeneeee 
problem of coupled conduction and free convection for a 
neatede triangular fin rotating together with surrounding 
air has been thoroughly studied. Second, by using this 
periodic heat transfer problem as a vehicle, a rational 
approach has been developed for solving this type of 
periodic problems characterized by finite fluctuating 
amplitudes in the solutions. 

AS sa result of rotation, three body forces are 
involved in this periodic heat transfer problem. The 
centrifugal force provides a persistent drive to the fluid 
flow while the gravitational force brings the periodicity 
dnto’ the problem. ~The Coriolis force Asean tay have any 
Significant effect on either the temperature or velocity 
SOmItLONS 7 1t can however contribute to the pressure gra- 
dient in the lateral direction of the boundary layer. The 
ratio, of thescentrifugal force to tneygrayvrcational force 
is the Froude number -(Fr) for thissrotating system. 

This Froude number was found to be the most important 


parameter in the present study. Formulation of this periodic 
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heat transfer problem revealed an unsteady-state, two- 
dimensionad;emathematic modeljoftesix control parameters: 

Ver Ge Kyo Ky, Fr, and Pr. By means of a Lefevre-type 

Emans formation nwethe effect OfPraon the solutions; was 
virtually eliminated. Among the remaining five parameters, 
the Froude number (Fr) would EP Re el tie fluctuating 
amplitude of the periodic solutions; the larger the Froude 
number the smaller the fluctuation. For an extremely large 
Froude number, a quasi-steady-state system containing only 
parameters y and C was obtained. 

The quasi-steady-state solutions of the simplified 
time-independent equations (3.7), (3.8), and (3.13) were 
studied for different values of parameters yand C. The 
resulbtsmag recduwellowath ‘Enose. of pEevilouseworkml2/7 2 sie 
In general, the parameter C had a significant effect on all 
the solutions while the parameter y could hardly influence 
the temperature field. | 


Two empirical equations in the forms of 


=O So 0G 


= a(C)y (4.5) 


10 
| 


and 


J 
G 
| 


= OSA ae (c)y]Rat/* (4.9) 


were obtained for the total heat transfer rate (aim) and the 
Nusselt number (Nia) at the fin root under this quasi-steady- 


state condition. Values of the functions a(C), Fy (Cc) and 
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F4 (C) were tabulated for various levels of parameter C. A 


least-squares fit of Oy iin ta simpler form of 


s a 


Of. Shah (@)y.4 

was also conducted and the function a'(C) was approximated 
by a fourth-order polynomial of the parameter C; the maximum 
percentage error of 10% due to this See fit was considered 
acceptable from a practical engineer's viewpoint. The 
tabulated values of Fy were also .of practical,use to.a.design 
engineer in estimating the periodic heat transfer rate from 
a heated rotating fin system (appendix E). 

Perturbation solutions were used to study the periodic 
problem for small fluctuations,and the investigation was con- 
1! Ko, and 


Fr. The parameter C had a significant effect on both the 


ducted for different values of parameters y, C, K 


velocity and the temperature solutions; whereas, the parameter 
y presented a greater effect on the velocity than on the 
temperature. The Seeioebe Gem seeuiles also revealed that in 
general, the periodic heat transfer problem of a rotating fin 
was nearly symmetric about the two vertical fin solutions 
(straight down and straight up). Namely, solutions at the 

two horizontal fin positions were normally identical. Only 
when K 7710 and C>>1l, did these two solutions have appreciable 
difference in their velocity profiles... The..fourth parameter 
Ky, however, had no bearing whatsoever on the boundary layer 
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the periodic phenomena were concerned, periodicity appeared 
Clearly in the velocity solutions when Fr=10 and became 
vanishingly small when Fr=100; no significant fluctuation 
was ever present in the temperature solutions. 

The perturbation results for a representative 


Sx cuationy Ofpe—0n0157 ve0e1y OK =100, were compared 


yqt-o, and K, 
with the quasi-steady-state solutions for the same values 
of y and C. The total heat transfer rate of the perturbation 
solutionss for Fr=10> matched closely with both the quasi- 
steady-state result and the mean value of the periodic total 
heat transfer rate corresponding to a smaller Froude number. 
The perturbation solutions also showed that the fluctuating 
amplitude of the total heat transfer rate stayed small when 
Fr=100 and increased to approximately one per cent of the 
quasi-steady-state result when Fr=10. Inspection of both the 
Nusselt number and the fin effectiveness revealed that their 
Brie Guat one were smaller than that of the total heat 
transfer rate. The average Nusselt number for any Fr value 
coincided with the quasi-steady-state result, but the average 
fin effectiveness seemed to decrease with Fr. 

A far more prominent difference between the periodic 
total heat transfer rate, the Nusselt number, and the fin 
effectiveness was in the phase angle. For the aforementioned 


sample problem, the maximum and the minimum total heat transfer 


rates occured at the 0° and ceo. position angles of the 
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rotating fin respectively; whereas, the maximum and the 
minimum Nusselt numbers at the fin root appeared to have a 
time lag equivalent to 10° angle behind these two extreme 
GOca. Neat, transter rates. On the other Nand, the maximum 
and the minimum fin effectivenesses took place approximately 
30° before the rotating fin reached these two vertical posi- 
tions (i.e. 0° and 180° respectively). 

The numerical solutions acquired from the direct 
integration of the time-dependent, governing equations 
(3.17)-(3.19) were used to study the periodic heat transfer 
problem for a small Froude number. Owing to the long and 
expensive numerical calculation of the coupled conduction- 
convection heat transfer problem at each time level, the 
study was mostly centred on the correlation between the Froude 
number and the periodicity. For the representative example 


Smee 04, YOU. L, hat, and K,=100, Dtawas, LOouUnds that 


He 
although the decrease in this most Significant parameter 

(Fr) did introduce some degree of periodic fluctuation to the 
temperature solutions, the oscillation was in general small. 
Tr contrast, ctoOetils, che variation of the Froude number 
appeared to have much influence on both the velocity profile 
and the heat transfer rate. When Fr=10, the directly inte- 
Grated velocity solutions OScil lated with respect to) time 


(or position angle of the fin) in the same manner as the 


perturbation result except that the Latter indicated a 3% 
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higher maximum velocity. Basically, this periodic variation 
in velocity profile remained symmetrical about the two 
vertical position angles (0° and 180°) of the fin until the 
wProude) number idecreasedto-ordénm of ones= int thatscase’,- the 
whole velocity Field not only fluctuated more intensively 
but also became unsymmetrical about the two vertical fin 
positions. For example, the velocity profile for 90° was 
different from that for 270°. 

The total heat transfer rate, the Nusselt number, 
and the fin effectiveness of this periodic heat transfer 
problem for a small Froude number were also studied. As 
predicted by the perturbation result, the maximum and the 
minimum total heat transfer rates obtained from the direct 
integration of the governing equations did occur at the chy 
and 180° position angles respectively, while the corresponding 
extreme values of Nusselt number at the fin root showed a 
certain delay in occurrence. In both solutions, there was 
arelis0r phase angle between two consecutive extreme values. 
Moreover, both the total heat transfer rate and the Nusselt 
number increased their fluctuating amplitudes but reduced 
their mean values when the Froude number decreased. The fin 
effectiveness calculated from the directly integrated solu- 
tion, however, revealed a periodic phenomenon with entirely 
different phase shifts. As the Froude number was reduced to 


order of one, the occurrence of the minimum fin effectiveness 


35 


was delayed but that of the maximum value was advanced, 

1.e€. a shorter build-up anda longer fall-off times in each 
periodic cycle. Notwithstanding all these periodic features, 
none of the fin effectiveness, (total heat transfer rate, and 
Nusselt number ever yielded a mean value deviating by more 
than 3% from the quasi-steady-state solutions. 

Applicability of the above solutions is basically 
limited by the necessary condition imposed on the Rayleigh 
number for a stable laminar free convection flow. If the 
HL COL Ra<10° for the laminar free convection flow in a 
constant body force field is also used for this periodic 
heat transfer problem and air is assumed to be the fluid, 
the obtained solutions can only be applied to the rotating 
system where centrifugal acceleration is not greater than 
50g, such as in a re-entry environment or a low speed centri- 
fugal machine. These results may be applicable to a high 
speed turbine if the actual upper limit on the Rayleigh 
number for a stable laminar free convection flow in the 
periodically modulated body force field is greater than boce 
This critical Rayleigh number can be determined by the sta- 
bility study, for which the results obtained in this study 
provide the basic solutions. 


on a rotating fin 


Based on the assumption of Ra<l0 
presents better cooling effect than a stationary fin does. 


According to equation (B10), the mean heat ‘transfer rate 
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POL -agrotating fin 1s higher than wthateoc. thesstationary 
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fin can be improved to approximately 2.66. times that of:a 
conventional stationary fin if it is economically and 
technically allowed to rotate the entire system to an 
acceleration of 509. 

Evidence has verified the reliability of the 
techniques and numerical programs utilized in this work. 
The close agreement between the present quasi-steady-state 
solutions and the results of previous work [27, 28] indicates 
that both the Blasius-Howarth type transformation and the 
numerical program developed for this quasi-steady-state 
problem are reliable. Further agreement between the quasi- 
steady-state solutions and the results obtained by directly 
integrating the time-dependent equations (3.17)-(3.19) for 
Fr=10°, confirms that the simplification of these time- 
dependent, periodic equations into the quasi-steady-state 
Cquiationse(s =. /)., (3.8) ,pand 3) for an extremely large 
Froude number is proper. That is, the normalization procedure 
adopted during the derivation of the governing equations is 
Sound: sn fact ~sthatwal) ienewCalcusdcederuncterTonsmandaticrr 
derivatives are of order one is another evidence of the pro- 
priety of this normalization procedure. The agreement be- 


tween the quasi-steady-state solutions and the time-dependent 
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results for an extremely large Froude number also assures 
the dependability of the numerical scheme in solving the 
periodic problem. This assurance is furthermore substan- 
tiated by the fact that these directly integrated solutions 
match closely with their corresponding perturbation approxi- 
mations when the Froude number is much greater than one. 

Two iterative processes were used in the numerical 
scheme of this study. All the quasi-steady-state solutions, 
the perturbation results, and the time-dependent solutions 
at each time step were obtained in a similar manner by 
alternatively solving the conduction and free convection 
problems. This iterative process converged when the tem- 
perature solutions of the conduction and free convection 
problems matched at the solid-fluid interface. On the 
Other hand, the periodic solutions were sought through the 
timewise iterations, which continued until all the velocity, 
temperature, and their derivatives at each point throughout 
the two-dimensional (X-Y) domain converged. 

The parameter C had a significant effect on the 
solution convergence rate of the iterative process between 
the conduction and free convection problems. The smaller 
the parameter C the faster the convergence. 

The solution convergence rate of the timewise itera- 
tive process depended strongly on the judicious selection 


of initial conditions and time level. The periodic solutions 
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for a problem under very small fluctuation was expected to 

differ only slightly from the quasi-steady-state solutions; 

however, when these quasi-steady-state solutions were used 

Bomtnerstancingssolutionseiat ¢=17/ 4) forest r-00F andthe 

Periodic squatrvons,A 3 Li)-$A(3.19) weresuniegratedsdarectly, 

solution divergence was found after numerous timewise 

iterations. In fact, neither could the rapid convergence 

be produced by using the perturbatvon results as the starting 

solutions for the periodic problem under moderate fluctuation. 

It was then imperative to develop a method of assuring accu-= 

rate approximations to the missing initial conditions at a 

proper time level. 

A rational approach, developed in the present study 
to produce rapid timewise convergence by providing good 
Starting solutions at the proper initial time level of the 
period, is described below: 

Step (1) The quasi-steady-state equations (3.7), (3.8) and 
(3.13) are solved. 

Step (2) The obtained quasi-steady-state solutions are used 
as the zeroth-order perturbation approximations; 
the.perturbationwresuits Pieand)P2,,usay,,£on) tne 
small and moderate fluctuations are generated 
successively. 

Step (3) Based on the perturbation result Pl, a proper initial 


time level is selected; the direct integraiton of the 


Step (4) 


Step (5) 


Lame-dependent ‘equations a(S be (3219). vsitthen ini — 
tiated after all the time-derivatives neve been 
approximated by Pl. 

Integrations in the succeeding time levels can be 
Cazrledvoutimaccordingly by) replacing ‘thes time= 
derivatives in the equations with finite-difference 
approximations; the periodic solutions Sl for this 
small fluctuation problem are obtained through 
timewise iterations. 


Once the solutions Pl, P2, and Sl are obtained, the 


proper initial time level and the associated initial 


conditions for the moderate fluctuation problem of 
interest can be determined by the rules given in 
SeCH1LOns3e2..2: throughutherwaltustratronjot some 
typical examples. Experience has shown that if the 
initial conditions and time level are selected in 
this manner, convergence OLmtnensobubionsias 


generally achieved within three timewise iterations. 


The above approach was proved effective and reliable 


im handling?) this) diffusion-type, periodic problems: under 


moderate fluctuations. The: total computing time required by 


using this approach to solve the periodic heat transfer 


problem for Fr=2 was approximately 75 minutes on an IBM 3607 Ge 


computer. 


For comparison, a numerical experiment was con- 


ducted by approximating the initial conditions with either 
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une quasios teady—-state solutions softnStepe (bitonethe pertur-— 
bation result of Step. (2), and integrating the time-dependent 
eguationsi (3.1 7)3( 3219) (for: Pr=27— Thestorumers approximation 
led to solution divergence while the latter produced time- 
wise solution convergence after 200 minutes of calculation. 
It was obvious that the approach developed in this study 
could save a great amount of computing time: 
Leesnouldsbe;notedsthat ithe vlomgecomputatron ofy 75 
minutes by using the developed approach is mainly due to the 
complexity of the coupled conduction and free convection 
heat transfer model at each time step; the solution conver- 
gence rate of three timewise iterations is in fact very fast. 
As far as the application is concerned, this approach 
can provide a general guideline for solving other periodic 
problems which occur in the’ physical world but remain analy- 
tically untouched because of their complex nonlinear pheno- 
mena aggravated by the finite fluctuation in the solutions. 
A typical example is the huff-and-puff steam flood process 
widely used by oil industry to recover very viscous oil from 
an underground oil reservoir. A certain amount of steam is 
injected into the reservoir to increase the formation tem- 
perature near the wellbore and lower the oil viscosity; the 
wellwis then shutéiin forma icriticalspersodto® time toyketetie 
reservoir fluid flow into the wellbore; steam injection is 


re-started after this reservoir fluid has been produced from 


my) 


MEERA tat 


mieesame well ye iThe erficiency iorethus stecnnaiquemis very 
Sensitiveitothe amountiof steampingected and» the! critical 
shut-in time. The physical models developed in the labora- 
tory have limited application value because’ of the radical 
changes in reservoir characteristics from well to well; the 
numerical models can easily accommodate any reservoir des- 
cription but suffer a major setback due to the great diffi- 
culty in solving the diffusion-type, nonlinear, periodic 
equations for this intermittent steam injection and oil 
production scheme. One major oil company has spent multi- 
miliion dollars on conducting an;80-well pilot project in 
the Cold Lake area to determine the optimum steam injection 
rate and the critical well shut-in time. Apparently, a 
great amount of money and human efforts will be saved if 
the periodic solutions for this huff-and-puff process can 
be obtained by using the approach developed in this study. 
5.2 RECOMMENDATIONS FOR FUTURE WORK 

The stability study of this periodic heat transfer 
problem is recommended. As shown in section 4.4, the 
applicability of the obtained solutions depends mainly on 


the critical Rayleigh numer for having a stable laminar 


free convectionet low. @ Pbfithescritical@Raylergh number deter— 


minedmoy  theystabidityestudysolthismperiodicihneatatranster 


problem is im thejorderyot 1078 the obtained solutions may 


bes applicablievto! theshighyspeed turbines. (inprecent years, 


fig 


RUusSiam SCientists have carried out) extensivesstability 
studies [16-24] on periodic free convection problems ina 
closed space. The periodicity enters the problems of their 
studies either by assuming a periodic temperature solutions 
or by vibrating the entire system in a gravitational force 
field. In most studies, high frequency modulation has been 
found to have a stabilizing effect on the ‘solutions. Should 
the same conclusion be deduced from the stability study of 
the present periodic heat transfer problem, the applicability 
of the solutions presented in this work would be extended 
due to a greater critical Rayleigh number. 

Experimental work to verify the theoretical result 
of this study is also recommended. The convection heat 
transfer phenomenon may be observed by means of a Mach Zehnder 
interferometer, and the temperature of the solid fin can be 
measured by thermocouples. In order to obtain representable 
interferograms, the relative distance and alignment between 
the fin and the optical plates of the interferometer must be 
exact to within 0.001 inch.. Therefore in designing the rota- 
ting system, the possible mechanical vibration should be 
taken into account. Moreover, the FM transmitter can be 
used to transmit the electrical signals between the rotating 
system and the stationary instruments in the laboratory. 

It is also worthwhile to examine the possibility of 


the flow reversal when L-y<a> eveboundary (layer tequationm a(3...3) 
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indicates that/if the icoefficient (1-yx) +> cos (277) of 
Chesbucyancyslorce term is positive-all) the time or, more 
specifically eave l-y>a> at the fin root, the fluid always 
tllows “in (the vradial inward direction from the tip to the 
BOOoEODNthe heatedvrotating fin. ~ For l-y<a>, the flow is 
sculLiIMinethemradval inwardldirection wien tthesrotatingnfin 
Zoepnuojeceing Iyertically tdownward  (7—0mand [eos (21t)=))). 
However, the situation may change when the fin rotates to 
the vertical upward position (t=1/2 and cos(27t)=-1); the 
buoyancy force tends to drive the fluid in the radial out- 
ward direction wherever L-yx<e> - Due to the inertia of 


the fluid, the flow may not reverse itself when 1-Yx is 
Bs 


Slightly less than Pr On the other hand if the reversal 
of the flow does take place, it would most likely start from 
the fin root because of the minimum value of l-yx (i.e. 1-y¥); 
that is, a new leading edge of a boundary layer, where the 
‘fluid-Llows radially outward, may EOuME at the pti mn loot 
whenever the rotating fin approaches the vertical upward 
position. Hor; a given geometric ratio y, there will be a 
eGritical valuervor Fr for the excitation of this reverse-flow 
phenomenon. \This critical Froude number is certainly worth 
exploring. 

A further investigation is to determine whether the 


reverse flow can really develop so that the flow will oscil- 


Vate valong the sfin’ with each» revolution, 0G to the ‘contrary; 
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the reversal of the flow will immediately break away the 
boundary layer. 

ineclosimng7 eit is Strongly ereconmencded tO apply “the 
developed approach to other periodic problems having non- 
Synchronous solutions of finite fluctuating amplitudes. 
The immediate, feasible applications are to solve those 
periodic problems presented in references [1-9] for moderate 


and large fluctuations. 
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APPENDIX A 


DERIVATION OF EQUATIONS 


Under the assumptions given in section 2.1, the 


governing equations of this rotating fin problem are: 


+ Up, + ees + pv) =—uOe CAS) ) 


OREO TS oe VUy) er bed, WU (Uyy + Uyy) 
(A.2) 


- p[ (R-X)w* + 2A0V + gcos(wt)], 


e(V, + UV. + VV.) 


x y -P. + H (Vy + V 


) 
i es (A.3) 
+ p[Yw* + 2AwU - Agsin(wt)], 


Re: 2 2 
pc, (Ty + UT, + vTy) a K, (Ty + Toy) + 2u (Uy + Vy) 


, 2 (A.4) 
+ HW (UY, + Vy.) P 


where subscripts t, X, and Y denote the partial differen- 

tiations with respect to these variables, and A=-1 or +1 

for the trailing and leading sides of the fin, respectively. 
In laminar-free-convection problems of incompressible 

flow, the temperature difference is usually regarded as the 

only significant cause of the variation in density and, also, 

this variation is assumed to be small. In other words, if 

8 is the coefficient of volume expansion and subscript "*" 


indicates a convenient steady state of reference, the density 
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of the-fluid at temperature T is taken as 
poe sil Ce Gea (Ae 5) 


where 
Bit =T a <<a (ARG) 


is the necessary condition. 

When the problem is steady, a conventional assumption 
is often applied; namely, the density p is regarded as a 
constant in all terms of the governing equations except 
those terms associated with body forces. This is because 
the body forces usually have larger magnitudes than the 
others (such as inertia force) in free-convection problems, 
the product of the small density change and these body 
force terms can not, then, be ignored in spite of the 
Sonditvons(A.o). Lniract, athis. produce as tie primary 
factor driving the free-convection phenomena. However, 
the feasibility of applying this Oberbeck-Boussinesg appro- 
ximation to an unsteady free-convection problem needs some 
further consideration. 

When the unsteadiness is due to the loss of 
equilibrium from a steady state such as in transient 
problems, the transient terms in the coverning equations 
can be small in comparison with the body force terms. 

The application of the Oberbeck-Boussinesgq approximation 
is then essentially the same as in the steady problems. 


In the present study, the unsteadiness is caused by the 
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periodical yavary ing wellectvon the —<GravutattonaleLorce. 
If the rotating speed of the system is relatively high 
but not high enough to abate the gravitational force 
effect, the transient terms may not be small and, on the 
contrary, may even be of the same order of magnitude as 
body force terms. Hence, it is possible that changes in 
fluid density can produce a significant effect in 
periodic convection problems through not only the body 
force but the transient terms. Therefore, substitution 


of approximation (A.5) into equations (A.1)-(A.4) gives, 


respectively: 
a eek ee UL IV) le ea ea (A.7) 
oes g (T-T_) (U, + UY, + VUy) = ras + yyy + Uyy) 

=—{l SA Wie te) (RES) ae 7 2have are geos(at) |, (A.8) 


[1 - g(t-T)\v, + Uy + vy)== P 


co 


+ v Vyy + V 


X X¥ WY 


+= eht HT) DY econ aur Scere ely (A395) 


a e(T-T ty UT. VEY |= Ker) 


».4 yan 


2v 2 2 Vig 2 
a aa lee) Mae =a Cou ed) Te (A.10) 


For convenience, the temperature of the fluid 
Minfinitely far from. thewtin ‘can pe aciiosen to beethe 
reference temperature of the system; in addition, it is 


also desirable to set 
Pos Pe ot P (Che MS ts) 


Le 


ee 
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where P is the pressure departure from the reference 


PESSSUrS sdiStriputlion Py vsatasiving 


oP 

5x7 = -o [ (R-X)w* + gcos(wt)] CA) 
and oe 

oS pf ¥w* - AgSin(wt)]. (ASS) 


PRmOCUCTEWOLCo yer (iS the pEessure distr loutionewhen the 
Peters orctvelcw Gt. AN aunTtOmmLenpetatunem ln .e.uned tier 
Heating nor cooling at the root of the fin) jand 1s rotating 
oem CONS Cantmangular velocity Wrabourmauruxed snori zontal 
axis. 

With the help of equations (A.11)-(A.13) and the 
condition (A.6), the governing equations (A.7)-(A.10) are 


Simplified to be 


Uy + Vy = B(T, + UT, + VTy), Cres) 
Ue UUs Se VU = lp + v(U te) oN OV 
c x iy. 0. —x xx YY 


co 


+ B[ (R-X)w* + gcos(wt)](T-T_), (A.15) 


ye il 
Ug yy nye tvWVyy + Vyy) + 200 
- B[Yw* - Agsin(wt)](T-T,), (A.16) 
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These are the boundary layer equations used in the present 
Study. It is worth mentioning ithat owing to condition 
(A.6), the buoyancy effect of the Coriolis force is 
suppressed in the equations of motion by comparing it 
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APPENDIX SB 


NORMALIZATION 


The method for separation of equations, described 
in section 2.4, can be used to split the governing equations 


(A.14)-(A.17) into a steady state set: 


Ss Sv. s,s s,s 
ea UO ja Oy), (Bia) 
Ss SS... =1 2s Ss s 2 s bes S 
nis Weve enon lic te yay) 2rAWV~ + Bw (R-X)O, 
(B22) 
s_.S SEAS ye a RS Ss Ss Sea oS 
Be Vo Rineg Poe AAU Sana el meme e GY Bue y Oa; 
(B53) 
s,s S'S s s 2vV S.2 S$. 9 
US, t VO, = kK, (0,5 + py) + Es CUS ney) 
v Ss sve 
+ ee (Uy + Vy) (B.4) 


and a periodic set: 
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Splitting the boundary conditions in the same 
Manner, one can first analyze the steady problem separately; 
and then, based on the steady results, study the periodic 
part of the problem. 
Bees STEADY EQUATIONS 

Since the problem approaches quasi-steady state 
when the rotating speed of the system becomes very high, it is 
reasonable to expect that the steady solutions are those 
of references [27,28]. Accordingly, the same normalization 
procedure can be applied to the present problem. 

Upon using the following normalized variables [274 
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equations (B.1)-(B.4) can be normalized: 


co a = pO (unO> + ve") (B. 10) 
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(Bales) 
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When air is the fluid, Prandtl number is, of order 
one and the conditions that Rane BOL, and Os are all 
much less than unity, all the right hand sides of equations 
(B.10)-(B.13) are negligible. The following boundary layer 


equations are, therefore, derived: 


uy inv =" 0, (B.14) 
Ss 
= (cua vous) = He + (1-yx) 6%, (B.15) 
p> = One (B 216) 
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u'o? + 0° 6° = O° (Br, 
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As an example, for a 6-inch aluminum fin, 12 inches 
from the center of rotation, rotating synchronously with 
Pies samo Went. aire 300 rpms wt Oe ern te is 10°F, one will have 
80.*0.019, Os=2.2x10*, Ramen Oh ote nae Eeee Raa PS 0MOs . 

Equations (B.14)-(B.17) imply that in the boundary 
layers: 

(1) the viscous force and the buoyancy force due to the 
centrifugal acceleration are of the same order of 
magnitude; 

(2) the inertia force has a comparable magnitude to the 
above two forces if Prandtl number is of order one; 

(3) the inertia force is negligible if Prandtl number is 
much greater than one; 

(4) the Coriolis force induces the lateral pressure variation; 

(5) conduction and advection are equally important; and 

(6) disregarding the pressure distribution, the problem is 
symmetric about the X-axis. 

In a word, the nature of this quasi-steady problem is 

governed by the simultaneous equations (B.14), (B.15), and 

(B17), ana the pressure distribution can be jcalculated 

separately, if desired, after the velocity and temperature 

distributions have been obtained. 
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x = 0: ue = 0, g= = 0 
y Orie enue OL = 0, 1 eee ee = 9° (x) (B18) 
y = ©: a=! 0; g° = 0 


s : , : } 
where 00 &) is the steady temperature distribution along 


ehevtin tsuriace. 
B.Z2 PERIODIC EQUATIONS 

In addition to the normalized variables (B.9) for 
steady solutions, the following variables can be introduced 


into the unsteady equations (B.5)-(B.8) of this periodic 


problem: 
p p 
U V Xx ny: 
a — Thal ’ ye — <7 , x — an ’ Be coon val , 
Cc Cc G Cc 
(B.19) 
p p 
ci Pp Q 
T= rete r pP = Di i ane oP= on ; 
C Cc Cc 


where subscript c indicates the undetermined characteristic 
values of the corresponding variables. The normalized 


equations, then, become 
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In order to simplify these unsteady equations 
through normalization,, one Nas to findythe: characteristic 
values of the periodic variables (B.19). The fin length 


is a convenient choice for Xos Thus 


Y Seat (B.24) 
G 


implies that 
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When»the angular velocity 0 issdecreasing from its 
very large value for the quasi-steady case, the periodic 
phenomenon gradually develops in the system. It is quite 
clear that this periodic character comes from the buoyancy 
force effect when gravity introduces its periodic influence, 
gcos (wt Tt) in the X-directional equation of motion, into 
the quasi-steady temperature field o>. In other words, 
Eersoperlodicldriving force, Bgo°cos (wt Tt), is Warge= in 
the vicinity of the solid fin because of higher temperature 
or) decreases gradually across the thermal boundary layer, 
and eventually dies out with o° outside this steady thermal 
layer. Hence the Y-directional characteristic length of 


the steady problem can be chosen as Yor ies 


Y= L/Rat/* (B.25) 


and thus 
A Sa 
A consequent question of interest is how the periodic 

buoyancy force, ¢° cos (wt Tt), gives rise Ea ae velocity 
fluctuations. . Equation (B.21) shows that the fluctuation 

in uP may be,sinduced, throughy inertia, viscous, orseven 
transient terms. Consider the solid-fluid interface where 
the no-slip condition eliminates not only the transient 


and inertia terms but also the viscous term uP 


of equation 
xx 


(B.21). One can presume that the driving buoyance o° cos (wt.t) 
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is balanced by the viscous force Oye inecacevyaecinity «of 
the interface and set 


a a 
BO g = VU /Y ; 


which yields 


There are two possible ways the periodic velocity 
uP can introduce the fluctuation into the temperature field; 
one through the continuity equation (B.20) and the other 
through the energy equation (B.22). In the continuity 
equation, ai. the oe term is he sean Desens” of the 
Pluctuationwand, renee possesses the highest order of 
Magnitude, it must be balanced by the oF, age or vo oP 
terms. The fact that the last two terms are of the same 
Magnitude leads to the necessary comparison of the 
magnitudes of the first two terms. To find their ratio; 
it he aera one should first determine the characteristic 
time, tor of this problem. 

It is evident that the period of the gravitational 
effect, gcos(wt), which generates the whole periodic 


phenomena, is a reasonalbe time scale, to: Therefore, 
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indicates that the aforementioned ratio becomes 
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if the ambient *Liutd asaireac room, temperature and 


oe is 10°F, the order of this ratio approaches O yn? SO LW, 


Henee an equatzon, (B20), both te and a” oF should have the 


Same order of magnitude if the temperature fluctuation was 


produced through this equation. By setting 
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Subs tir bution “Olyenis oo value into the transient 
term oF in the energy equation (B.23) yields 
ae o 
a oP 2h re: 
Obviously, this should never be greater than the order of 
magnitude of BOE Which) Us tne source», Of the Liluceuation 


in this energy equation. However, the ratio of the magnitude 


of uP oe to oF reveals va contradteli on: 
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when air “isthe fluid ate room temperature, 

This conflict indicates that the temperature fluc- 
tuation does not come directly from the continuity equation. 
Nevertheless, one does find that the transient term oF is 
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It is then reasonalbe by equating ue and v® to get 
= ve 
1 
K-Ra ie 
Ma = hk Reale (B.28) 


With the aid of characteristic values (B.24)-(B.28), 


the energy equation (B.23) can be rewritten as 


| oP + uF oP + vegP + Sie (Oe Corer vi 02) 
DRat’ cit x y Ex G x Y 
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where parameter D is the dimensionless period for the fluid 
Elow, and is equal ‘to 2mK -/wL*. As mentioned before, if 
the temperature at the root of the fin is 10°F higher than 


miyatwOfagcSamolent aim; 1.e. O =10°F, then 


1 saa 
prai/2 27 LBO,. 
=o [y}/2] 


Since the parameter y is always less than unity, it is apparent 


that if A=DRai/? 
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— = < 
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In general, if eee Paine Os are much less than order 
one, all the viscous dissipation terms in energy equation 
become negligible in comparison with the driving terms, 
uP ® 
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of fact, the whole energy equation may be normalized by 


setting 


therefore the characteristic value of oo is around 


(BieZo)) 


Substituting characteristic values (B.24)-(B.29) 


into the continuity equation (B.20), one obtains 
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No contradiction takes place by choosing the definition 
e229) for ee 
In order to determine Par one again introduces the 
obtained characteristic values into, and then normalizes, 
the Y-directional equation of motion. Equation (B.22) becomes 
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Comparison of the magnitudes of parameters 
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will indicatenwhethermthe Coriolisvor ithe gravitational 


force causes the pressure fluctuation. The ratio of 


Sy to S.5 is defined by 
S = nih E bee ee 
S 21 Pee / 
ioo7Oltt, setting 
Po = pak pwRal/*/20Fr (B. 30) 


into equation (B.22a) gives 


pe = 2duP + [=]\6°sin (207) (B.22b) 


encdmrr 550 (lL) > setting 


Po = 0 ,B0_.gL/Ra‘/* (B. 30a) 
into equation (B.22a) gives 
Py = \o%sin(2mt) + [S]2AuP ; (B.22c) 


for. S=O[1l}. . Both definitions) (8.30) end, (8B. 30a): are iproper. 
in fact, nowmeatter whichrdccmntction/ oF Po is chosen, 


the X-directional momentum equation (B.21) can be normalized 
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= uP + 6 °cos(2mt) + (1-yx) oP + [==] oP cos (2 nt) 


~ Opal’ * Ek Ra’! *] Pr(Be2i 2) 
With the help of the example illustrated in section 
B.l, it is easy to show that the velocity and temperature 


fluctuations are governed by 


uP + vP = 0, (B.20b) 
x y 
eas S15 ps Sp Ps Lae ae P,,P 
Pra! 4 + u uy + u uy + Vv uy + Vv Bh + [a] (u uy + Vv a 
= uP + bd ®cos (2rt) + (l-yx)oP + [+ j6Pcos (207), (Be 215) 
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subject to the boundary conditions 


ae (Oe SRS OE oP = 0 
Veta 0 (reo) uP = 0, vP = 0, oP = oF (x,1) 
SE ue = 0, oP = 0 (2531) 


Pp 


The pressure filuctuationsp can be obrained separately from 


either equations (8. 22b) or (B.22¢) aderending= on, che magnitude 


of parameter S. 
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From this normalization, equation (2.8) can be 


re-written as: 


ep yi) eed oy a) er [aclu (x,y,t) + 

Vil x ee) acts) + [acl vP (x,y,t), (BuS2,) 
an’ ¢tyy,t) = oS Gey) + [el 6? ry, 7), 
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v= uae 
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are the normalized solutions of the complete problem. 


Accordingly,. the normalized forms of equations (A.14)-(A.17) 


are Wwe towle= 07 (B, 34) 
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ole and Je are not normalized guantities if the magnitude of 


W/Prods notyofyorcer onessi nce ft romecnce equalities §(B 52) 


and 


Pp Ag ee D 
AG >. Fr ee 


Cakes 
ee Ala a 


Hj 


\ ‘ te 
ean! Pa. Ae 
j 
ea sae 
; a art, 4 


m 


i 
ms 


168 


where quantities GH and oF instead of ae and he , are 
actually normalized. This fact would not be revealed if 
equations (B.34)-(B.37) were obtained from the direct 
normalization of equations (A.14)-(A.17). In other words, 
when the Froude number Fr in equation (B.35) decreases 
with the angular velocity w of the system from oe high 
magnitude, this equation (B.35) can not provide sufficient 
information about the point where the transient term 
Seanrts introducing a comparatively sSsigniticant effect on 
Bhe problem. |/From’ equation {By21b)) tt is very clear that 
the parameter A is the determining factor of this. 

If the temperature at the fin root is 10°F above 
the ambient temperature and the Prandtl number of the 
fluid is of order one, the magnitude 


ieee as 1/2 


Let z 1/2 
x= tome =e Oey sacs (Bs) 


indicates that the transient effect is negligible when the 
geometric ratio yY approaches the limiting case of y=0 
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APPENDIX C 


LOCAL SIMILARITY TRANSFORMATION 


In order to solve three simultaneous boundary layer 
eecat1 Ons. (B254))) 1 41(B. o5)) anda (hao). amstreame function 
Wwex, vy) derind by 


. ov me 
Esom oy and Bas ee y 


Uselatrogdiced. BSincerthe fcontinuity equation (B. 34) is 
Satisfied completely, this boundary layer problem is then 


governed by only two simultaneous equations: 
L 
P + ~ + Pr j(1l- ales cos (211) 
Ue ay ele ial ty es [ ye lee 9 


ie ee = 0 Coors) 


ee 
byt Vyby ~ Wyo, - igld, = 0 Ges 


The associated coundary conditions are: 


aC; by ae Lely, d= 0 
y = 0 Cie xe Oye: p= 0, vy = 0; od = a (x) ’ 
Vie tees Vy =a 07, d= 0. 


The follow-up ettoOrtei ssi ae tandingpasurcabie trans— 
FOEMALION (CO; 
(1) lessen the dependence of the solutions in the streamwise 
direction and thus decrease the error introduced by any 
finite difference formula on approximating streamwise 


derivatives of the governing equations; 
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(2) alleviate the possible singularity problem at the 
leading edge; that is, 


Lim. (x) +2 1imd(07y)o=.0 
x20 yO 


(3) reduce the importance of the Prandtl number in the 
resulting problem and thus remove the need to determine 
Sseparace Solutions .tor a/warge number of Pr. 

A transformation is assumed to take the following form: 

& = L(x), re ete sore) ANIM orgy EG) 
(Cer } 
®(€,n) = $(x,y), and E(e ea Vey) / Gt Pay Mix), 

where (&,n) and £(&,n) are two dependent variables in the 

new coordinate system, and correspond to original temperature 

function $(x,y) and stream function (x,y) respectively. 

Se un knowl LunctLons, oi (x)),, M() 5 1 yy ee Cer) mana G (Pr), 

are disclosed after the substitution of this transformation 

miconeduations, (CL) “and “(C.2).8 That as, 

(per tcun*34)£,,,, + 3(PrF°GMN*IDI)£,, + (PFET yyy) £,| 
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Wisehiithespurposeconsitructingmesghacis (fon) local 
Similarity", terms containing only the n-derivatives are 
kept on the left-hand-sides of the equations. In addition, 
terms for the same physical effect are square-bracketed 
ito ONeHGLroOup. -yhoOr imstance, the tirst Ehree groups on 
the left-hand-side of equation (C.4) are, ina sequential 
order, related to viscous, inertia and buoyancy forces 
while the first and the second square-bracketed groups 
on the right-hand-side correspond to the variation of the 
inertia force along §-direction and the transient effect. 
Similarly, the first and the second groups on the left- 
hnad-side of equation (C.5) represent advection and 
conduction heat transfer while the first and the second 
right-hand-side groups indicate the variation of conduction 
in the €-direction and the transient effect, respectively. 

It is assumed at this stage and is later substan- 
tiated by the normalized equations that the highest order 
derivative term in each square-bracketed group should never 
be suppressed by other terms in the same group. Hence, 
onlyJthe coefficient of Highest order term nia group is 
employed below to determine the aforesaid six unknown 
functions whenever the group is considered to be physically 
Ssigniticant. 

Since advection and conduction are equally important 


in the thermal) boundary Layer, tte reasonable to equate 
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their parenthesized coefficients in equation (C.5). That 
oe! 212 = FGNM_I 

yy; x y 
or by rearrangement, 


N (x) G (Pr) 
MY Gat) ty) TE (Py) val 


where Ci is an undefined constant because the left-hand- 
Side of the above equality is a function of x and y while 
thesright-hand=side is/a ‘function of Prvonly. Similarly, 
N(x) _ ee 5 
Maina Sippel ie koe: 


that is, the function I(y) takes the form 


ie 


G 
I(y) = ome! +C 


For convenience and without loss of generality, one can 


goeune C_=C,=1l and C.=0 and obtain 


Lies 2 3 

B( Pa). = G{Pr) (C6) 
La) ee a, | (era7) 
Rs em ee tc) | (@S)) 


iImpordersto find B(Pr) -OrsG (Pr) cons: deration 
should be given to the momentum equation (C.4) where the 
buoyancy force term, %, is always dominant. However, the 
inertia LoOrceacerm, eae and the viscous: force term, Ponte, 
vary inversely in significance as the Prandtl number 
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forces play) important roles’ an balancing’ the buoyancy force. 
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CPA ee Pate a. 
[F “GMN I pr=o[1] = Ouel| 
yields 
3 4 o 
MM LF Ipreof1] = Ol1) (C399) 
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When Pr approaches infinity, the inertia force is negligibly 


small but the viscous effect remains significant. That is, 


MMe dim [Piao (Cell) 
r7o 
and 
344 Pa 
MM iim Pr = 0: (Cou) 
Pr-o© 


On the other hand, the viscous force can be dominated by 


Enesanertia, Lorcevac Very ysmal., Pr. manus, 


MM?lim [F*] = 0 (Gus) 
“Pr7*0 
and 
3 F 
MM, lim lF| = Onl). (C22) 
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are completely satisfied. This can be fulfilled by setting 


ee as 4 ie 
MM, a Aly and FE 7 Pr =i). 
Hence, 
SCD 3/ 4 
Mix) = (5 x + C,) 
and 
re Pr 1/4 
aleoet rome(Ci 8) 
4 -1/4 
N (x) See) = (> x + Cy) ; 


Since: in transformation +(C23.), 
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the possible singularity at the leading edge of the fin, 


ane’ lim , (x) > 0 
x70 


oie lim ¢(0,y) = 0, 


gral 


can be removed by setting C,=0 to get 


et fae 0) y. 


In other words, the singular point (0,0) in the x-y plane 
is "stretched" into a’ line in the,new'é—-n plane so that 
the discontinuity of temperature at this ssingulan point 
is ‘smoothed out. along the line. 

In} order) to define ‘al proper tz unction, for Lix)\vone 


can substitute the following functions: 


4 3/4 -1/4 
L(y) =; M(x) = (= x) ; N(x) = (= x) 


and 


Ud 


F(Pr) = G(Pr) = (GEE) 
antLoreduations (C.4) Vand, (CG. S)mtouget 
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(Cos) 
To complete the transformation from the x-y domain 

into the €-n domain, all the square-bracketed terms in the 
above two equations should be either a function of & or 
constant. Assuming 

ee) ay 
yields 

—E = L(x) = yx 
Waerelthe! aroitrary, ConstantsoL integrations set etoube 
zero for the sake of convenience. To sum up, the governing 
equations (Cl); and (C..2))are transformed) into 
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by the local similarity transformation: 
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The associated boundary conditions for the above equations 


and 


Ee ty) 


are 


oe Wa EO ems RON )G 1 S1(57/0) 
(@x2 1) 
QN = ©; ie L = 0, $= 0. 
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APPENDIX D 


PERTURBATION EXPANSION 


The purpose of adopting this technique in the 


PEeesentu Study 1s moe only (tovexpand, thersoliitions, £(677),0) 


ny Oe, Nyt) soo equations Ase! /)—(3..9),.1n a perturbation 


series of 1/Fr, but also to remove the time function from 


the perturbed equations so that less computation time is 


required to solve the remaining two-dimensional problem. 


To accomplish this, one seeks solutions in the form: 


£(E,n,t) = £°(E,n)+ SAlE(E,n) cos(2mt) + F(E,n) sin (207) ] 


+ Sec aa =. £(E,n) cos (47T) 
Fr? 


+ ee un) sin(4nt)] + aa 
Fr? 


and 


(Deas) 


B(Emn,t) = OF (En) + FlS(E,n) cos(2nt) + $(E,n) sin(2nt)] 


a 1 5(E sn) f 6(E,n) cos (47T) 
EG 


+ $(€,n) sin(4nt)] + OL] 
Bre 


where f° and ae the quasi-steady solutions when Fr= 


satisfy the zeroth order equations: 
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The first order perturbation solutions £(§,n), 


£(E,n), O(E, 


Hospeand 6(E,n) etc. can be obtained from the 


following simultaneous equations: 
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with the corresponding boundary conditions 
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yields solutions.f and’?, etc: 


The second order perturbation solutions f, 
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and 6 etc. should satisfy the following simultaneous 


equations: 
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The associated boundary conditions of the above equations 
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Tt is clear that the time variable has been removed 
from the problem. As a result, difficulty in’ predicting 


accurate initial conditions (at every single point on the 
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Sono ULLdce)mcOnasSure TLast convergence sor the speriodic 
solutions is avoided and excessive computation caused 

by time-wise iteration is effectively eliminated. In 
Spite of the fact that this approximation technique will 
inevitably become tedious and less reliable when much 
higher-order perturbation solutions are required (for 
the small Froude number cases) its lower-order solutions 
alone still provide some extremely valuable information 


and hints in solving the problem for small Froude number 


tsecelon 392.22). 
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APPENDIX E 


A SIMPLE TIME DEPENDENT SOLUTION 


When the thermal inertia effect on the phase 
behavior of solutions is neglected, the Nusselt number (Nu) 
for laminar free convection around the heated EOtatingg@e rn, 
depicted in figure (2.1), can be estimated by a simple 


engineering approach. At first, the following relation 


is assumed: 


1/4 


Nu « Raa'ce (Ba) 


where the effective Rayleigh number, Ra segs is defined by 


3 
COR rere 


Races = cay (Bi2)) 


If the streamwise acceleration at a point on the fin surface 
L/n feet (n21) from the fin tip equals the effective 


acceleration Tore for the entire fin, then 


Sog¢ = (R-L/n)w* + geos (wt). (E.3) 


By using y=L/R and Fr=Rw*, equation (E.3) becomes 
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Tore (1-y/n)Rw? + gcos (wt) (E.4) 
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Furthermore, combination of equations (E.2) and (E.5) yields 


g{(l-y/n)Fr + cos (wt)] gO L* 


cee Jane VK ee 
2 
one: 
Races = [(l-y/n)Fr + cos (wt)] Ra_, (S37!) 
where 8g0 _L* 
Ra. = ara ° (E38 


From equation (E.7), relation (E.1) can be rewritten as 
Nu « [(l-y/n)Fr + cos (ut) ] a eee (E.9) 
Or more specifically, 
Nu = m ((l-y/n) Fr + cos (wt)]1/*Ral/* (B70) 


where m iS an empirical constant. 

The value of m can be estimated by considering the 
Nusselt number at the fin root for the quasi-steady case 
ep in 


equation (E.2) is replaced by the streamwise acceleration 


of the problem. When the effective acceleration Se 


atethe rin rootu n=.) ppeduatiton) (hr. 20s Cannbe used ato 


calculate the associated local Rayleigh number; that is, 


B [(R-L)w* + gcos (wt)] 01° 


Ra, (t) = Th T Ome RINERO WEEE ee Cee ly) 


Consequently, the local Nusselt number at the fin root can 
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DewoxpDLesseo Dye SubStLtutingun=laintowequation. Wwm.10)e) namely 
Nu, (t). = m[(l-y)Fr + cos (wt)] 1/*Ral/* : (Be 12) 


When the problem is in the quasi-steady state, the Froude 
number (Kr) Gn ‘equation, (6.12) asso jJarge that the: cos (at) 


term can be neglected. Therefore, equation (E.12) becomes 


NO = m [(1-y) Fr] ?/*Ra?/* “ may)? “(erRa,)?/* 
hee aly Rates « (E.13) 
where BRw*0_L? 
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By using the series expansion: 


(1-y) 1/4 = 
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kK 
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equation (E.13) is rewritten as 
N& = m(1-2 Ra?’ ‘ | (fan ale 


Where ‘constant () .LS-used)to account sfor, Che error or 
approximation. 
Lt acpworth) poanting, outethat®’equationgs(h. 14) eas 


in the identical form as the empirical equation (4.9): 


Na = (Eons Fiy)Ral/*, | (4.9) 
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which has obtained independently in section al ey, by 


applying a least-squares curve fitting’ technique to the 


quasi-steady-state solutions. 


(4.9) and (2.14) suggests that 


m= Fy CEA Se) 


Values of Fy are tabulated in Table (4.4) for 
various levels of parameter C, which represents the degree 
Of coupling between the conduction and convection heat 


transfer problems. 


Once the value of m is known, the time dependent 


Nusselt number for the rotating fin problem can be determined 
by equation (E.10) provided that the thermal inertia effect 


on the total heat transfer rate is negligibly small. Fora 


sample problem of C=0.01, y=0.1, and Fr=2, Table (4.4) and 
equation (E.15) yield 
me= 0239-0 


and equation (E.10) becomes 
Nu = 0.39[2(1-0.1/n) + cos (ot)] var ane (E.16) 


If the streamwise acceleration at the middle point (n=2) 

of the fin surface is chosen as the effective acceleration 
Sore for the entire fin, the time dependent Nusselt number 
e 


is then expressed as 


Nu = 0.39[1.9 + cos (wt)] */ "Ral? * : (Eee 


Comparison between equations 
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The reliability of the above approach can be checked 
as follow. The maximum and minimum values of the Nusselt 
humber, corresponding to t=0 and t=1/w respectively in 


equation, (E.wL7) ~.are 


Z 
Gc 
II 


0.5089 pale 


and 


Z 
G 
II 


0.3799 Baas : 


the ratio between these two extreme values is 


= 1.340 (E.18) 


The dimensionless total heat transfer rate based on the 
theoretical solutions for the same example is shown in 
figure (4.27), where the additional parameters Ky and K, 
regulating the thermal inertia effect are set to 1 and 100 
respectively. The fact that the Maximum and minimum values 
in figure (4.27) also appear at t=0 (a=0°) and t=m/w(a=180°), 
indicates a negligible thermal inertia effect on the total 
heat transfer rate; that is, solutions are in phase with 

the periodic acceleration gcos (wt) for this sample problem. 
Despite the thermal inertia effect, the ratio of the maximum 
to the minimum total heat transfer rates presented in figure 
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Now, because the total heat transfer rate is proportional 
to the Nusselt number, the ratio obtained in (E.19) on the 
basis of the theoretical solutions can be compared to the 
met OLestimaceduin (BE. le) Mbyausing joisepracticalapproach? 
a difference of 2% is found between these two ratios. 

This observation suggests that equation (E.10) aie Table 
(4.4) can be used as the practical means to estimate the 


total heat transfer rate for a heated rotating fin problem. 
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